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COMPRESSIBILITY DATA FOR HELIUM AT 0° C AND PRESSURES 
TO 800 ATMOSPHERES 


by 


Ted C. peigee =! He uss faitorel! and Robert E. perteate! 


ABSTRACT 

Twenty-two compressibility runs were made with helium at GAG 
and pressures to 800 atmospheres by using a Burnett~-type apparatus. 

The purpose of making multiple runs at one temperature was to 
statistically determine an adequate functional form for representing 
the data, and to then present precise new data. 

The work is part of a broad program for the experimental determ- 
ination of data for the evaluations of the thermodynamic properties 
of helium and helium-containing mixtures. 

Data from the 22 runs were treated simultaneously and as sep- 
arate runs while different functional forms and weighting factors 


were used. A summary of the results is included. 


Research chemist, Helium Research Center, Bureau of Mines, 
Amarillo, Tex. 
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A fourth degree polynomial equation in the pressure gives a 
statistically satisfactory representation of the compressibility fac- 
tor when the data are treated simultaneously and a weighting factor 
of prt is used. 

The second virial coefficient of helium at 0° C from this inves-~- 
eAgationsisi lie9970e.0.015 cm? mole +. The stated uncertainty is the 
random error standard deviation. Compressibility factors determined 
by this investigation are precise to better than 0.01 percent and are 


believed to be accurate to better than 0.02 percent. 


INTRODUCTION 

Accurate experimental compressibility data for helium and helium- 
containing mixtures are of importance to the Bureau of Mines helium 
conservation program. These data are needed for calculations of the 
amount of helium transported through the various pipelines, of the 
amount stored in the conservation fields, and of the amount ultimately 
delivered to the user. Thermodynamic properties, computed from the 
compressibility data, can be used to design improved process and sep- 
aration cycles eae the extraction, concentration, and purification 
of helium. 

Extensive and accurate compressibility data over wide temperature 
and pressure ranges are required for the development of an accurate 
equation of state. These data will contribute toward planning and co- 
ordinating the utilization of the nation's helium resources, and will 


aid in the future development of new uses for helium. 
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The Burnett compressibility method (12)3/ was selected for the 


uy 


=" Underlined numbers in parentheses refer to items in the list of 


references preceding the appendices. 


experimental measurements. This method does not require volume cal- 
ibrations, nor is it necessary to know the number of moles of sample. 
nee precise experimental observations can be made with this method. 
One objection to the Burnett method is that it does not yield an 
experimental compressibility factor directly. A functional form must 
be assumed before the compressibility factor can be derived from the 
experimental observations, or the data must be reduced graphically. 
We feel that a generally acceptable treatment of Burnett data, 
free from esoteric manipulations, will be required to obtain agree- 
ment among investigators consistent with the experimental precision. 
Criteria must be developed for the selection of a functional form to 


represent the data and for the choice of a weighting factor. 


it is toward these coals that »sthis.work is directed. 
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EXPERIMENTAL APPARATUS AND EXPERIMENTAL PROCEDURE 
The apparatus used in this investigation was a conventional 
Burnett compressibility apparatus (12). A Burnett apparatus basically 
consists of two volumes designated as V, and Vz, equipment to measure 
the pressure of a gas sample confined in volume V, or (V, + Vg), and 
equipment to control and measure the temperature of V, and V2. 
An unscaled diagram of the gas pressure system of our apparatus 


is designated figure 1. Figure 2 is a drawing of the oil pressure 


FIGURE 1. - Gas Pressure System of the High-Pressure Compressibility 


Apparatus. 


FIGURE 2. - Oil Pressure System of the High-Pressure Compressibility 


Apparatus. 


System. Figure 3 shows the constant temperature bath along with 


FIGURE 3. - Constant Temperature Bath System and Some of the Major 


Components of the High-Pressure Compressibility Apparatus. 


some major apparatus components. 

Some of the equipment incorporated into the present apparatus 
has been described in the literature (14, 29, 36, 38, 47). A detailed 
description of the component parts of our apparatus is contained in 


appendix A. 
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PyGUR@ec —Oilehressure System of the High Pressure 


Compressibility Apparatus. 
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FIGURE 3.—Constant Temperature Bath System and Some of the Major 


Components ot the High-Pressure Compressibility Apparatus. 
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The Burnett compressibility method has been described in the 
literature (12). 

Some variation in experimental techniques is possible as long as 
the data treatment and experimental procedures are consistent. Con- 
siderable attention to detail and careful experimental work are re- 
quired to obtain maximum precision and accuracy from the equipment; 
therefore, a detailed description of the experimental techniques used 
in this investigation is contained in appendix B. 

CALCULATION OF ABSOLUTE PRESSURES FROM THE 
EXPERIMENTAL OBSERVATIONS 

A computer program was written by the Data Processing Branch to 
calculate absolute pressures from the experimental observations. 

The first step of the calculations was to evaluate the barometric 
pressure from the observed scale reading and temperature. National 
Bureau of Standards Monograph 8 (11) was used as a reference for re- 


duction of barometer readings to standard conditions. 


We take 
R R I+st C 
< = < ( b) ey 7 ’ (1) 
760 760 \l+mt,/ g, 760 
where 
R, ° 
760. barometric pressure, standard atmospheres 


R, = barometric scale reading, mm Hg 


v 


‘ - ty 
_— 


4 
xt 
is 
sii 


s = coefficient of linear expansion of barometer scale, 
PGsatl0re. Ci... (khy..page 28) 
m = mean cubical coefficient of thermal expansion of 


mercury 
mxl0° = 18144.01 + 0.7016t, + 0.0028625t2 + 0.000002617¢3 , 
if eels epage 4) 
ty = barometer and/or room temperature, °C 
& = acceleration due to gravity at the barometer 
location, 979.4091 cm sec *, (18) 
8, = standard acceleration due to gravity, 980.665 cm sec # 


C, = calibration correction provided by the National Bureau 


of Standards for serial No. 13346 barometer, +0.06 mm Hg 
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Next, a piston gage pressure was computed from the experimental 


observations. National Bureau of Standards Monograph 65 (17) was 


used as a reference for reduction of piston gage data. Equation 2, 


_ Ma (1 - Pa/Py) & /Se 
Fs ~ A, (+bP, )L1te(t,-25)] 7S > ae 


was used to calculate gage pressures for the data recorded in this 
report. Calibration constants and data for piston gage, serial No. 
9274 were provided by Ruska Instrument Corporation. 

In equation 2, 


P, = gage pressure, psig 


xe el 
ee om 

7 q Y 

7 ss i ii 
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ct 
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M, = mass of the piston gage eae! , lbs 


The calibrated piston gage weights were reported by Ruska 
Instrument Corporation as apparent mass in air against brass 
standards; therefore, the standard brass density of 8.4 grams 


cm * was used in equation 2 rather than the actual density of 


the weights, (17, page 3). 


0, = density of air, grams cm ° 

Py = density of brass standards, 8.4 grams cm ® 

g = acceleration due to gravity at the piston gage location, 
979.4091 cm sec *, (18) | 

Z. = Standard acceleration due to gravity, 980.665 cm sec * 


A, = piston gage No. 9274 effective area at 25° C and 
atmospheric pressure, 0.0260416 in® 
b = piston gage No. 9274 effective area pressure distortion 
coefficient, -3.5xl0 ° psi? 
c = piston gage No. 9274 effective area temperature 
distortion coefficient, 1.7xl10 ° °c? 
t, = piston gage temperature, °C 
Cy = differential pressure cell correction, 1.85xl0 * P,. 
Barieau (1) derived equations 3 and 4 for the calculation of the 


density of moist air in grams liter ?: 


D 
1 
{ 
wyelq ais to tena 
ne Sees eS Om 
| rs ris 
1 TU ‘ 
; 
f 
2) 
fa 
ee 
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> ; 
; 
! 
j 
- 
' ce Pa 
1 > y t 
_ we 

, erry | 

, i ; 7 

> n = ys 

re : 7 P. 


pr #152932 = — (1-0.37807y) 2 , (3) 
<2 = 1-0.000602 (= 3 1) + (0.00254y + 0.0758y") se 
760 
-(0.0000105 + 0.0000131ly + 0.00131y”)t, a : (4) 
where 
Ps = density of moist air 


Te@raveolute temperature, 2/3.15° + t,, °K, 
T, = absolute temperature at the ice point, 273.15° K, 
and 


y = mole fraction of water vapor. 


Equation 5 was assumed to be valid for the air density calculations 


of this work: 


, MLE ee 
where 
P,; = vapor pressure of Liquid water, atm, and 
H = relative humidity expressed as a decimal fraction. 


Osborne, Stimson, and Ginnings (42) presented the following 


equation for the calculation of the vapor pressure of water. 


3 3 : 
log; oPs = -3.142305(A07 : Eta?) + 8.2 ae 
-0.0024804(373.16-T) , (6) 


where T = (273.16 + t,), °K. 
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The density of air was calculated by using equations 3, 4, 5, and 
6, the barometer or room temperature, barometric pressure, and rela- 
tive humidity. 

Piston gage pressures in psig were calculated from equation 2 by 
using an iterative procedure which continued until the change in the 
computed gage pressure was less than 0.0001 psig. Calibration data 
for piston gage No. 9274 along with the experimental observations 
were used as input data for the calculations. 

Defined values used by Ruska Instrument Corporation in the cali- 


bration of piston gage No. 9274 were: 


Vw ele, 2.54 cm exactly 


isi: 


453.5924 grams. 
An atmosphere in terms of the constants supplied by Ruska Instrument 
Corporation is: 


+ 1.013250-10°:(2.54)* _ . 
1 atm = 453.5924 980.665 14.69594780 psi 


Ruska piston gage pressure in psig was divided by 14.69594780 to ob- 
tain gage pressure in atmospheres. Absolute pressure in atmospheres 
was obtained by adding the barometric pressure in atmospheres to the 


piston gage pressure in atmospheres. 


EXPERIMENTAL RESULTS 
An experiment was run to determine if any impurities were intro- 


duced into the test gas while filling the compressibility volumes to 


ee 
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the initial high-pressure. An eee e eno cylinder of grade A helium 


a 


Lys 


=’ All composition analyses were performed by the Branch of 


Laboratory Services, Bureau of Mines, Helium Research Center. 


was attached to the inlet of the gas compressor. The gas system was 
evacuated and purged with test gas. Then the pressure containers 
were filled to an initial pressure of about 700 atmospheres. Fill- 
ing of the pressure containers was carried out exactly as it would 
be done for a compressibility run. Gas was vented from the pressure 
containers, and a sample of the test gas was trapped for analysis. 
The compressibility volumes were filled three times to pressures of 
about 700 atmospheres, and three samples were trapped for analysis. 
A sample taken directly from the supply cylinder and the three vent 
samples were analyzed by high-pressure mass spectrometry (48). The 
results of these analyses are recorded in table 1 of this report. 
Examination of the data of table 1 indicates a slight increase in the 
nitrogen and oxygen content of the vent samples over that of the cyl- 
inder sample; however, the increase in total sample impurities is 
considered negligible and probably occurred during the sampling pro~ 
cedure rather than during the bomb filling procedure. 

An experiment was conducted to determine if any errors were Jie 
troduced into the compressibility pressure measurements due to cdif- 


ferences in the density of oil in the line to the differential 
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TABLE 1. - Analysis of vent samples arts per million mole 
fraction impurities in helium 


Sample direct from 


cyiindér No. 25598... 0.2 
Vent sample No. l...... az 
Vent sample No. 2...... < 

4 


Vent sample No. 3...... 
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pressure cell. Part of the vertical portion of the oil line was at 
ambient temperature while another vertical portion of the oil Line 
was at the bath temperature. This situation existed because the oil 
line entered the constant temperature bath over the edge of the bath 
rather than directly through the bath wall. 

A heat exchanger was placed around the vertical part of the oil 
line located outside of the constant temperature bath. Bath fluid 
could be circulated through the heat exchanger so pressure measure- 
ments could be made with the outside vertical part of the oil line at 
room temperature or at approximately the bath temperature. 

Measurements were made with both sides of the differential pres- 
sure cell open to atmospheric pressure, with and without cooling in the 
heat exchanger. A measurable difference in pressure due to oil den- 
sity changes was observed. This effect was completely compensated 
for at atmospheric pressure by setting the meter of the differential 
pressure cell on zero at the start of a Poaeeesalbility run while 
atmospheric pressure was applied to both sides of the differential 
pressure cell diaphragm. 

Next, pressure measurements were made at about 365 atmospheres 
and at about 700 atmospheres with and without cooling in the heat ex~ 
changer. There were no measurable differences in the observed pres-~- 
sures when cooling was applied to the oil Line or when the line was 


at ambient temperature. 
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We conclude that any errors in our pressure measurements due to 
density gradients in vertically displaced tubing lines in inhomogeneous 
temperature environments are completely negligible. 

A total of 22 compressibility runs were made with helium at O5eC! 
Three helium cylinders were used to supply sample gas for the runs. 
Analyses of impurities of each cylinder are recorded in table 2. The 
particular cylinder used for each run is listed in table 3 of this 
report. 

A platinum resistance thermometer was mounted in the compressibil- 
ity bath and was used to measure the temperature of, the bath fluid (9). 
The average bath temperature for a complete run was assumed to be rep- 
resentative of the actual sample temperature. The average measured 
bath temperature and the standard error of the average measured bath 
temperature are recorded in table 4 for each of the’ 22 rungaec idem 
peratures are in terms of the International Practical Temperature 
Scale (IPTS). 

Experimental gage pressures, barometric pressures, and experi- 


mental absolute pressures are listed in table 5 for each run. 


TREATMENT OF THE EXPERIMENTAL DATA 
The data treatments are based on the mathematical formulas 
developed by Barieau and Dalton (2) which are specific applications 
of their general nonlinear least-squares techniques (3). Two computer 
programs, written by Douglas and Moffitt (19) and by Grout (22), were 


used in the data treatments. 
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TABLE 2. - Helium cylinder analysis arts per million mole fraction 


impurities in helium 


Cylinder 
number 


H- 135022 
AEC107250 
AEC 103941 


TABLE 3. - Helium cylinders used to supply sample gas for the 
22 compressibility runs 


Cylinder No. 


Cylinder No. 


“ H-135022 


0-1 AEC 107250 

0-2 H-135022 AEC107250 

0-3 H-135022 AEC107250 
HE-0-4 H-135022 AEC107250 
HE=0-5 195022 AEC107250 
HE-0-6 H- 135022 AEC107250 
HE<0-7 H-135022 AEC107250 
HE-0-8 H- 135022 AEC107250 
HE-0-9 AEC 107250 AEC103941 
HE-O~10 AEC107250 AECLO3941 
HE-O-11 


AEC107250 AEC103941 
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TABLE 4. - Compressibility bath temperatures 


HE-0-1 .006 +0. HE-0- 12 +0.003 +0. 
HE-0-2 +.002 +.001 HE-0-13 +.002 +.001 
HE-0-3 -.001 +.001 HE-0-14 +.003 +.001 
HE-0-4 -.001 +.000 HE-0-15 +.002 +.000 
HE-0-5 +.001 +.001 HE-0- 16 +.003 +.001 
HE-0-6 +.003 +.001 HE-0-17 -.001 +.001 
HE-0-7 +.002 +.001 HE-0-18 +.002 +.001 
HE-0-8 +.004 +.001 HE-0-19 -.001 +.001 
HE-0-9 +.003 +.001 HE-0-20 .000 +.001 
HE-0- 10 +.002 +.001 HE-0-21 +.002 +.001 
HE-0-11 +.001 +.001 HE-0-22 .000 +.001 


TABLE 5. - Experimental pressures 


pressure ,atm D 


REAR S _Run No. HE-O-1 
0 .4120277069 0.8769311813517 
1 280 .3759539201 .8769454385195 
2 130 .8076096956 .8770905151851 
3 62.97670429211 . 8768062323970 
4 30 .61519494207 .8764649185141 
5 14 .78507020371 .8762032717097 
6 6 .946035181388 .8756657411545 
tee ees -039906680540 |  .8752078666919 
ae Sa et) wake ee Oe Run No. HE-0-2_ 
0) .5170315508 0.8762738697152 
1 296. 7309855857 . 8763450974200 
2 137 .9581306295 .8760239146634 
3 66. 32798378899 -8750997728367 
4 3224010065330 .8747869576745 
5 15 .58767440123 .8745395418001 
6 . 345869800917 .8738342433041 
i? TE a oe) 5240439740641. | .8735213919652 | __ 
a ee eevrrene sees Ron Noss hEeOs3o)4y) 
0) .4084415110 0.8739253759817 
1 300 .8341453471 -8738514027795 
Z 139 .7399746708 .8737916238975 
3 67.16019770479 . 8733621694687 
4 32 .64421980891 .8731005311728 
5 15.78791517340 .8726568518288 
6 -446612599290 .8719885684152 
EG ae -292375244210 | _.8716102816285 

iat. eS Jerenrio eu Run NocaHE-0-4 
.8641147696 


302 .4895877704 
140 .4576848864 
67.49191770212 
32 .80078704206 
15 .86074344101 
477834751585 
3 301370424995 


INAUFWHFRO| 


Run 
0. 


- 9955688057 
299 .5308764732 

- 1709653445 
66 .88916039177 
32 .50828178547 
15.71605782455 
-406375775771 
eOGS0GOLObbL« | |. 


NOUR WHE CO 


Barometric 


ressure ,atm 


8817856682931 
8818367398739 
8817345949371 
8813483717485 
8810642718613 
8812318696604 
8812686284498 


o. HE-0-5 ef) 
8839887329134 


- 8838435598045 
-8834511534051 
- 8830077009650 
-8819102020566 
-8812561701240 
- 8804488692424 
8796332827422 


Absol 
_pressure,atm _ 


. 2823668870 
L. 7079967498 
6137662947 
03355987426 
51732034008 
66057202523 
318601167705 
163985525839 


. 2889588883 
2528993587 
6847002107 
85351052451 
49165986059 
66127347542 
821700922543 
__3.915114547232 
3933054205 
. 6073306831 
8341545442 
20308356183 
. 11488761097 
46221394303 
219704044221 
113961132606 


7452014818 
. 3713734387 
3395216263 
37365229705 
68213541380 
74180771287 
359066621246 
182639053445, 


8795575386 
4147200330 
0544164979 
.77216809273 
39019198753 
59731399467 
286824645014 
147941902353 
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TABLE 5. - Experimental pressures~-Cont inued 


pressure ,atm 


Run 


1184744399 
4683901262 
0183182667 
29454910138 
.71106911175 
15..82199250306 

7 465224361390 

3 .302930725135__ 


IINUF WHR CO 


2180830037 

. 1394439427 

- 8803547994 

- 22943089199 

- 68057988169 
15 .80835835797 
458410384007 
. 300149583854 | 


INDUFP WHF OC] 


Run_ 


Absolute 
pressure ,atm 


Barometric 


“hu 1 0 eee 
.8702066886767 702 .9886811286 
.8700616745118 302 .3384518007 
.8700333791464 140 .8883516458 
.8691492883398 68. 16369838972 
.8687656722470 33 .57983478400 
. 8687887124937 16.69078121556 
-8682371602237 .333461521614 

___.8680691486822 | 4. 170999873817 

_No. HE-O-7 Remenin © Pee were 

0 .8687210533628 702 .0868040570 
- 8686824068501 302 .0081263496 
. 8685233086120 140. 7488781080 
.8684065799029 68 .09783747189 
.8674251969547 33 .54800507864 
8666259990989 16 .67498435707 
.8664156377437 .324826021751 
—WBODTATIS2Z67522 fh . 165896910606 
No. HE-0~-8 


0 7285681673 8827347709460 6113029383 

it 300 .5655371206 8829246892377 301.4484618098 

2 139 .6251276605 8831145738312 140 .5082422344 

3 67. 10290632809 8825912546956 67 .98549758279 

4 32 .610143324705 8826199551148 33 .49405320216 

5 15. 76649427462 8823009046303 16.64879517925 

6 7 .430836014840 8815383662152 312374381055 
boven 278450102214 | _.8813788877451_ 159828989959 _ 
Beasts Oe OO Se em 
0 683.7113986306 0.8633966156373 684 .5747952462 

1 294 ..5442035018 #663550/333729 295.407 7542352 

2 137 .0068278494 8636442251013 137.8704720745 

3] 65 .889699 73099 8632890950770 66.75298882607 

4 32 .0333363/603 , 8632236855030 32.89656006153 

5 15 .48782008142 8627658184852 16 .35058589990 

6 7 .300724831583 8622285236672 . 162953355250 
ead | 223765520253 | _.86 18734563292 1. .085638976582_ 
‘hey ee ee YT un Nos sHE-O0=10 ERIS ST TE OO 
0 .4006329458 0.8749632389877 707 .2755961848 

if 303 .0742241189 .8756968986971 303 .9499210176 

2 140 .7169389485 8757396096306 141.5926785581 

3 67 .61494253624 8759785524275 68.49092108867 

4 32 .86278893303 .8761890262546 39. 3097/95929 

5 15 .89149483270 -8760070366065 16.76750186931 

6 7495482000912 .8757453898021 .371227390714 
7 _3.312819082249 8754837429978 188302825247 
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TABLE 5. - Experimental pressures~-Continued 


Gage pressure ,atm 


-5241244213 
6190721062 
.0822349999 
6/ .J4/01oIoI2d 1 
32 .71699099093 
1>.619996/3512 
7 460373454598 


NOU WHF CO 


1774194280 
9914171748 
3770283460 
6699035375297 
32 .56209784277 
15 ..74527942310 
7 424333849166 
279510932827 


NaUPWNF OO] 


- 1568450183 
-4851924201 
- 2684818409 
67 .91060783425 
Lan OM AOE RS a te Vs 


15.97185770344 
7.541110926879 


INInUP WNP Ol 


- 6627960235 

wh POD LOS 

- 2935027453 
66 .95088667224 

32 .54100831655 
15 . 73488840020 
7-418015287165 


NSounPWN Fe O| 


40013 
»0579222501 

- 7030908260 
67.60547615335 

32 .85474805965 
15.88603113031 
7490382185862 


NOW WN Fe CO 


0.8794799381871 


.294951232367 | 
Run_No. 


Run No. 


» 342758393301 | 


-275455464607 | 


- 308926021058 | _ 


Barometric 
pressure ,atm 
HE-O-11 


-8796925959585 
-8795617789376 
-8793798502741 
- 8785723770434 
.8779182493966 
- 876609994 1030 


HE-O-12 


8753698827309 
8750713191034 
8750201394122 
8736750129967 
.8731005311728 
8725772545810 
8723668427877 
HE-0-13 
0 .8676085841898 
| 8671403682317 
8669672182456 
| 8668081781562 
8661026033712 
8654997101146 
8648454760937 
8648314140170 


Run No. HE-0-14 


. 8760808845412 


~8759785524275 
»8753386879669 
ehh 9 L272, 49.13 
8746988025848 
30/4582 1963979 
Run No. HE-0-15 


. 8818250796224 
.8817820635829 
.8816225369381 
. 8812299837589 
.8807576837478 
. 8803508223873 


8763995068130 


"S758277050837 || 


0.8761605363940 | 


Ler SHA | VANS M7 I is Wa 


0.8821074330062 


_ .8799582819748 | —- 


Absolute 
_ pressure ,atm_ 


4036043594 
4987647022 
9617967788 
68 . 19639374298 
3359556336798 
16.69791500511 
336983448702 
171350739181 


0532471331 
8667870575 
2520996651 
67 86537389239 
3343577285577 
1661837995427 
296911103747 
Rel LOTT SG Wise 


..0244536025 
hs ip aah yad Keio we 

. 1354490591 
68.77741601241 
33°..0/7 (13453219 
16.83735741355 
-405956402972 
207589807318 _ 


-5 389565599 
6756759550 
1695324647 

67 .82686522467 
3341634700452 

16 60997967299 
292714089749 
150037663005 _ 


2539614343 
9397473297 
5848728896 
68.48709869029 
33.73597804341 
16. 76678881405 
. 370733008250 
. 188884303033 
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TABLE 5. - Experimental pressures~-Continued 


Barometric 


Run No. HE-0-16 


0 8780447545 0.8743716205323 
1 302 .5013284971 8744598322460 
2 140 .4678649356 8741184531733 
3 67 .50087124635 8733418957295 
4 32 .80882900305 8724629627553 
5 15 .86828180075 8717661818688 
6 485697576303 8713082255012 
7 ~310740668733 _ 8705744244341 


_Run No. HE-0-18 


Run No, HE-0-19 


_Run No. HE-0~-20_ 


0 3298625462 0 .8756664095909 
1 295. 1451253806 8756976258424 
2 137. 2620746086 8760475950205 
3 65.99994771761 8754158428421 
4 32 .08007116705 . 8751910392938 
5 15 .50704145733 .874.7414747694 
6 .305086175129 | 8742037435632 
7 218774229971 | _.873 7969042076 


Absolute 
_pressure ,atm 


pressure ,atm _| 


0 2372753729 0.8750600359184 
1 300 .0107262375 8752220961934 
2 139 3833027406 8753102257581 
3 66.99246918307 8750201394122 
4 32.56161610313 8745054240785 
5 15. 74537885842 8736779874653 
6 423746309458 8731406014611 
7a 2781751 


8728364418387 | __ 


0.8726024839403 | 


0 1153690775 

1 302 .21899700i7 8726710398511 

2 140 .3451528107 8730863364027 

3 67 44352544444 8730067345339 

4 32..78196173269 8712265691693 

5 15 .85640853209 8697363621858 

6 480714618288 8694093249510 

7. .306957570342_ 8692926427963 __ | 


0.8758770341266 | — 


0 9865344622 

1 298. 7833283840 .8753618584328 

2 138 .8467393769 .8756748981965 

3 66.73951310222 -8759736808505 

4 32 .43839773162 .8757630779786 

5 15 .68246073581 -8762808882956 

6 - 390886624099 -8/591L09L91L181 
ee 3260606042031 |  .8762238875722 _ 


686. 2055289558 


705 . 7524163750 
303 .3757883294 
141.3419833888 
68.37421314208 
3368129196580 
16 . 74004798261 
8 .357005801804 
4181315093167 __ 


300 . 8859483337 
140 .2586129663 

67 .86748932249 

33 43612152721 

16 .61905684589 
8296886910919 
4.151011592503 _ 


704.9879715614. 


303 .0916680416 
141.2182391471 
68.31653217897 
33, 65318830186 
1672614489428 
8. 350123943239 
_4.176250213138 _ 


695.8624114964. 


299 .6586902425 
139 eI 224142751 
67 .61548678307 
33 ..31416080960 
16.55874162411 
8.266797543218 
_4 136829929603. 


296 .0208230064. 
139314661222036 
66.87536356045 
32 .95526220635 
16.38178293210 
8.179289918692 
_4 092571134179 _ 
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TABLE 5. - Experimental pressures~-Continued 


solute 


Ab 


‘Barometric 
Gage pressure ,atm 
42.000 eras | _ pressure ,atm 
_Run_No. HE-0-21 


0 479244748 0 .8746931169914 774 .3539378655 

1 327 . 9503923646 8746817427037 328. 8250741073 

2 151.4663134605 8743119936602 1523406254541 

3 72.62586671643 8737460381948 7349961275463 
4 35 .28527915768 8731289217587 36. 15840807944 
5 17 .08701089108 8725118450070 17 .95952273609 
6 091134431307 8721564377234 | 8 .963290869031. 
ie 612440477369. 8717498202373 _| .484190297606 


_Run No. HE-0-22 


0 . 3092872243 0.8766132786184 582. 1859005029 
1 255.2070832114 . 87649084 16644 256.0835740531 
2 119. 7052428099 8765361723173 1205817789822 
3 5775117285393 8766099941030 5862778284803 
4 28 .07558827206 8763198475297 28.95190811959 
5 13 .53087369724 8759643108626 14 ..40683800810 
6 320927090012 «| ~—.8758704223976 7.196797512410 
was -726188056005 | _. 8756457310696 | 601833787075 
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Data from the 22 runs were treated in different ways. First, 


each run was treated separately. Values for the parameters and 


compressibility factors were evaluated for each run when the data 


were fitted to the equations; 


2, 


and 


Zz, : 
1+BP,+ch =— £,NP, , (7) 
8 
: Zn Oe 
iter oednChsettOProcis- fyN Beso» (8) 
“0 
1+ BP, + CPE + DPR + EP, = 5 f,NP, (9) 


compressibility factor at P,, 

parameter evaluated from the experimental pressures, 
parameter evaluated from the experimental pressures, 
parameter evaluated from the experimental pressures, 
parameter evaluated from the experimental pressures, 
pressure after the r th expansion, 

pressure before the first expansion, 

compressibility factor at P,, 

R = expansion number, 


factor to correct for elastic pressure distortion 


isothermal zero pressure volume ratio. 
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Weighting factors of P,, l, pieek P,+, and P;* were used when 
the data were fitted to equations 7 through 9. A summary of the re- 
sults is presented in table 6. Details of the calculations are on 
file at the Bureau of Mines Helium Research Center. Data for run No. 
HE-O0-13 were not consistent with the data of the other runs and were 
not included in the calculations. 

We decided upon the following criteria for an acceptable data 
treatment: 

1. The solution should give pressure residuals that appear to 
be consistent with our assumption of a normally distributed random 
variable with a mean of zero. 

2. The solution should give values of W,(68P,)* which are the 
same for different expansion numbers except for statistical varia- 


tions, where 


= 
i} 


. weighting factor, 


§P standard error of a single pressure residual. 


P 


None of the solutions met these criteria when the data were 
treated as separate runs. 

Next a method was developed to treat all of the data simultane- 
ously. 


The function 


M 
Rie || dew, (aP,)? | (10) 
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Weighting 
factor 


FFP FP re 


.0005 36309 
000538871 
.000541274 


.000542158 
000546131 


.000533759 
.0005 34337 
000535614 
.000537852 


000525873 
.000526830 


.000528390 
. 000530567 


+0. 


000001334 
000001226 


.000001513 
.000001698 
.000002860 


. 000000336 
000000389 
.000000518 
. 000000930 


000000138 
.000000179 
.000000221 
000000314 


PF FP 


FFP FP re 


1347259 
5347385) ee 
254/633 =. 
9 SUIS IEE  a 


.345515 
13459866 = 
. 346442 
. 347075 


.347589 +0 
. 347832 
. 348033 
. 348100 
. 348353 


+ 


1+ 


each run separatel 


measurements on helium at 0° C, treating Pp y 


.000183 
.000179 
000211 
.000228 
000325 


Zy 
000072 
000086 
000112 
000174 


000044 
000056 
000064 
. 000099 


1.994197 +0.000052 
1.994121 +.000048 
1.994059 +.000053 
1.994038 +.000056 
1.993961 +.000076 


=12,.Po t,N Py 


5.36388 +0 
Blacteiles) as 
Dips tey2 es 
Sa 2) 2s 
By ova eh 28 


01341 
LOLZS3 
AQIS ZA. 
.01708 
.02879 


(Pep, te Ore + DP = 2,P, fp Ps 


1.994277 +0.000023 
1.994239 +.000028 
1.994171 +.000035 
1.994087 +.000049 


=rl-+ BP. HCP? = 
1.994812 +0.000019 
1.994659 +.000022 


4 


1.994459 +.000021 


1.994276 +.000034 


Hess) 240. 
Bin SEMIOY, Be 


5.35682 
By, SZ 
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lar 


Zoho EN By 

5.25919 +0. 
S266 lim 2 
5 284385 +, 
55306162 +. 


00337 
00390 
.00519 
. 00932 


00129 
00149 
00199 
00314 


11832) 
0 38 U 
AWE: 
1416 
4245 


.63860 + 
65164 
. 68357 
. 74848 
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- 46789 
47887 
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TABLE 6. - Summary of average quantities obtained from the nonlinear least-squares treatment of the experimental compressibility 


D:10°° ,atm ° | E-10!° ,atm * 
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is minimized in the simultaneous treatment of the data, 


where 
R’ = the sum of the weighted squares of the residuals of P,, 
M = number of sets of data, 


ny = number of expansions, 


i age 


where 
P, = observed pressure, and 
Py = calculated pressure. 


Each element of the coefficient matrix is calculated for each 
run, just as in the separate method (2), then the contributions from 
each run ate added together to obtain the coefficient matrix for the 
Simultaneous data treatment. Once this coefficient matrix has been 
obtained, then the rest of the calculation is the same as that used 
to treat the data as separate runs. 


Data from 21 ate were fitted to equations 7 through 9 with 


6/ 


= Data from run No. HE-0-13 were not included in the calculations. 


the method developed to treat data simultaneously. Weighting factors 
of 1) Poi Mees 8 P,*, and P;® were used. A summary of the results 
is included in table 7. Details of the calculations are on file at 


the Bureau of Mines Helium Research Center. 
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TABLE 7. - Summary of the results obtained from the nonlinear least-squares treatment of the experimental compressibility 
measurements on helium at On C, treating all runs simultaneously 


Weighting 
factor Z(1 atm Z(700 atm N B:10*,atm ~ C-10° atm = D:10°1,atm ° E-10!° ,atm™* 


= chee = 
7, = 1 + BP, + CP? + DPF EF, = Zaba) ENP 


1.000533574 +0.000000851 1.347212 +0.000154 -994285 +0.000046 | 5.33638 +0.00807 -6.4129 +0.2473 Ue 7 SOR Sa 2 =8:,33 
pi 1.000534673 +.000000701 1.347421 +.000108 994225 .000033 | 5.34740 +.00674 =6, 7415 2.2226 2 LOT = Sees =5,.83 21.93 
o™ 1.000535209 +.000000696 1.347517 +.000112 .994199 000030 |~5.35278 +.00672 -6.9099 +.2316 2,443 — bg37h mi) Lome Oo 
PS. 1.000535750 +.000000741 1.347610 +.000114 994174 +,000029 5.35820 +.00720 «7.0848 +.2592 2.701 £6425 -8.53 £2.44 
la 1.000537856 +.000001518 1.347944 +.000185 994092 +.,000041 Bg eyieys ssn AL) =i 09 Cae OL Sy AO koe elk IL IAG, -13.56 +7.09 


= Ce] = = 
Z. = 1+ BP, + Ge? + DP? = Z,P,,£,N PF, 


i 1.000532301 +0.000000288 1.346984 +0.000072 1.994349 +0.000025 . 32361 +0.00258 -6.0067 +0.0446 
1.000532779 +.000000272 1.347114 +.000062 1.994305 +.000020 NO26 39" eRO02 aI -6.0838 +.0500 
p;* 1.000533097 +.000000291 1.347195 2.000063 1.994279 +.000019 33159 eB, 00272 -6.1400 +.0580 
he 1.000533470 +.000000331 1.347284 +.000068 1.994253 +.000020 £33532 eee 00313 -6.2096 +.0710 
tie 1.000535361 +.000000751 1.347680 +.000124 1.994148 +.000030 £35427 +.00734 -6.5791 +.1941 


Z, = 1+ BP, + 
1.994946 +0. 


CPe=.2.,8,° E:N EF; 
000055 | 5.24580 +0 


1.000524535 +0.000000433 1.345098 +0..000137 -4.5120 +0.0240 


p 4 1.000526112 +.000000351 1.345707 .000108 1.994708 +.000040 5.26158 +.00302 -4,6130 +.0218 
ihe 1.000526792 +.000000323 345959 000098 1.994616 +.000034 52ZOSa9 -4.6588 +.0214 
7 1.000527421 +.000000309 1.346184 +.000091 1.994537 +.000030 5.27468 -4.7027 +.0219 
P;* 1.000529564 +.000000449 1.346900 +.000110 1.994321 +,000028 5.29613 +.00430 -4,8631 +.0413 
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The functional form of Z, selected as best representing the ex- 
perimental data, treating all runs simultaneously, is a fourth degree 
polynomial in P,. This particular functional form was the lowest de- 
gree polynomial for which there appeared to be random scattering of 
the calculated pressure residuals. Higher order polynomials in the 
pressure were investigated; however, no significant improvement in 
the randomness of the pressure residuals resulted from this work. 

The weighting factor was selected as follows. 

The product W, (6P,)* was examined for each fourth degree polyno- 
mial fit in the pressure to determine if W,(6P,)* was a constant 
independent of pressure. Standard errors of the single pressure resid- 
uals were fitted to the equation. 

W, (6P,)* = a + bP, Git) 
by using a least-squares method. Values of a, b, and b/a, along with 
their standard errors are listed in table 8. The calculated curves 
for 


W, (SP, )* 


b 
3 Sd ae, (12) 


are illustrated in figure 4. From table 8 and figure 4, we see that 


W, (6P, )* 
: a 


FIGURE 4. aseaerunction of P. . 


the most horizontal curve is that obtained when a weighting factor of 


a was used. 
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TABLE 8. - Least-squares values of a, b, and b/a calculated 
for the equation 


W, (8P, )* = a+ bP,, assuming 


ceoee Breet Ure + BPP + EP. 


5.764 +0.355 | 0.9064 40.6854 


pe 231 +.055 0585 +.0222 
Pia +. .003 +.038 .0014 +.0166 
ee +.399 -.034 +.031 -.0244 +.0190 
pe +.123 -.010 +.009 | -.0463 +.0342 


_6\d Sos .¢ .6 30 souley e9teupe-Jeaak = .o 2 


Wr =| 


SO 


iit oats 


1OO 


W,(6P 


a 


200 


35 


Wr = Pr (3 


3 
Wr =P, @ 


Zo) 300 


as a Function of P,. 
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Now we examine the fourth degree polynomial solution to determine 
if it appears to have normally distributed random pressure residuals 
with a mean of zero. 

Table 9 gives values of the pressure residuals as a function of 
the expansion number for a fourth degree polynomial representation of 
the compressibility factor where a weighting factor of pe has been 
used in the calculations. 

Table 10 lists thecalculated mean pressure residuals as a func- 
tion of the expansion number and the average pressure. 

The contents of table 10 indicate that the calculated mean pres- 
sure residuals differ from zero by less than two standard deviations, 
and only two of the seven differ by more than one standard deviation. 

We conclude that the nonlinear least-squares calculations give a 
solution consistent with our original criteria for an acceptable 
solution when: 

1. All runs are treated simultaneously, 

2. the compressibility factor is expressed as a fourth degree 

polynomial in the pressure, 

3. and a weighting factor of pv is used. 


Values of the parameters for the acceptable solution are: 


B = (5.35278 + 0.00672)-10 *, atm?; 
C = -(6.9099 + 0.2316):10 °, atm *; 
D = (2.4479 + 0.3709)-10 7+, atm °; 
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TABLE 9. — PRESSURE RESIDUALS, IN ATM, FOR A QUARTIC FIT IN P 
TO REPRESENT Zs TREATING ALL RUNS SIMULTANEOUSLY, 
WEIGHTING FACTOR P SUB R TO THE —-3/4 
R=1 


DEV. FROM AVG. 


RUN NO. (P,0BS.-P,CAL.) PRESSURE RESIDUAL 
ij So” deo 8 504355E-03 301 STEE—OS 
ne ee 2e6220E-03 2-9442E-03 
HE-O-3 teZ2385E-03 7 e5606E-03 
HE-O-4 -4.0043E-03 —SeDocLC—OS 
HE-0-5 —F~5102E-O3 -7.2481E-03 
HE-O-6 -2.8604E-03 “Levoore OS 
heE-o—7 Seto 1 EOS 5..6592E—-03 
HE-O-8 —2208165E-04 4.0463E-05 
HE-0~-9 “Lele -OS —2e030TE—-83 
HE-0-10 823244E-04 1.1545E-03 
HE-O—11 1.f422E-03 220643E-03 
He=-o-i2 3e.372TE-O3 3-6949E-03 
2BE-O-13 —5.23815E—-04 -2.1602E-04 
HE-O-14 5 -5341E-03 528563-E-03 
HE-@-15 —Letoroe—Os -9.1068E-04 
He O16 -222948E-03 i. 97275-05 
REO f 323049E-04 6.5261E-04 
HE-O-18 2e6247E—-04 5 -8460E-04 
nE-o-19 =o. t oor U5 Hie tt oe OS 
He 8-20 ~4.7740E-03 -4.4519E-03 
He-O-21 “Geb PDZE-OS -7.8907E—-04 
HE-@-—22 —T¥eFLOZE-O3 =f eSB LEAOS 


AVERAGE PRESSURE RESIDUAL -3.2212E-04 
STANDARD ERROR OF THE AVERAGE PRESSURE RESIDUAL Tn 3 SL SES 
STANDARD ERROR OF A SINGLE PRESSURE RESIDUAL 4e.2857E-03 


% DATA FOR THIS RUN WERE OMITTED FROM THE CALCULATIONS. 
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TABLE 9. — PRESSURE RESIDUALS, IN ATM, FOR A QUARTIC FIT IN P 
TO REPRESENT Zy TREATING ALL RUNS SIMULTANEOUSLY, 


WEIGHTING FACTOR P SUB R TO THE -3/4 
R=2 


DEV. FROM AVG. 


RUN NO. (P,OBS.—-P,CAL. ) PRESSURE RESIDUAL 
HE-O-1 4-8717E-03 4.6943E-03 
HE=O-—Z 328456E-03 3.6681E-03 
RhE-O-3 661297E-04 4.3548E-04 
HE-0O-4 —-1L.6425E-03 -1.8200E-03 
HE-O-5 2 «60,1 DE-O3 —20f TE6E-G3 
HE-0-6 —~Ze2ob9E-O3 =3..5089E-03 
HE—O—7 4.3843E-03 4.2068E-03 
HE-0~-8 4.6921E-03 4.5146E-03 
HE-O:-9 +344 50ZE—-O3 “Se SOE OS 
HE-@-10 Be2atSE-O2 220440E-03 
HE-O—-11 2.4180E-03 2e2405E-03 
HE-O0-12 321586E-03 2-e9811E-03 
*HE-O-13 -604667E-02 —-6 .4844E-02 
HE-O-14 4.6756E-03 4.4981E-03 
HE-6-245 =2.0612E=-63 =2e2381E-O3 
HE-@—16 ~3.5648E-04 he SOI TEES 
HE-O~-1T b.«Z571E-05 -1.6491E-04 
HE-0-18 4.9079E-04 3.1330E-04 
HE-0-19 -~4.5044E-03 -4.6819E-03 
HE-O-20 -4,.0421E-03 ~4.2196E-03 
HE-O~-Z1 -224658E-03 —226433E—-03 
HE~O—22 —Fa SODLE-6S ~3.4787E-03 


AVERAGE PRESSURE RESIDUAL 1.7749E-04 
STANDARD ERROR OF THE AVERAGE PRESSURE RESIDUAL 1 21064E-04 
STANDARD ERROR OF A SINGLE PRESSURE RESIDUAL S#ZOGsE-Gs 


* DATA FOR THIS RUN WERE OMITTED FROM THE CALCULATIONS. 
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TABLE 9. — PRESSURE RESIDUALS, IN ATM, FOR A QUARTIC FIT IN P 
TO REPRESENT Z, TREATING ALL RUNS SIMULTANEOUSLY, 
WEIGHTING FACTOR P SUB R TO THE -3/4 


R=3 
DEV. FROM AVG. 

RUN NO. (P,O0BS.-P,CAL.) PRESSURE RESIDUAL 
HE-O-2 3e-1747E-03 FA (oe pe od war | fe 
ne -O—2 2296 28C-65 2-6074E-03 
fe —o—3 2-681 7E-05 —-3.4864E-04 
HE—-O-4 —3.9924E—66 -3.1945E-04 
2B -1.4619E-03 4 -«8573E-03 
HE-0O-6 2e3657E-05 —-3.5180E-04 
HE—-O-/ 229341E-03 225586E-03 
Hi-8-8 3-1807E-03 228053E-63 
He—0—9 321092E-04 -6 .4544E-05 
HE-0-10 8.2330E-04 424783E-04 
HE-6—7 4 1.0911E-03 Te1570E-04 
HE—O—i2 220431E-03 1.6676E-03 
*+HE-O-13 226268E-03 Laz ot4E=-03 
HE-0-14 LeSSE9E-03 22457LE-03 
HE—O—15 —1.63536E-03 —h..fZ9LE-03 
HE-O-—46 726693E-06 =Je6 f TIE—S84 
HE—— 3 7 =J0 59 TSEL04 =t~eSl25E-0% 
HE—-O-—18 522692E-04 1.5145E-04 
Me—G-19 -328432E-03 -4.2186E-03 
HE-0-20 20659: 1LE—-93 —302346E-03 
ee! ome 1.2377E-04 —L0d LOBE—O4 
Me 0-22 —2+43200E-@3 —~2.6954E-03 


AVERAGE PRESSURE RESTDUAL Be 1546E-04 
STANDARD ERROR OF THE AVERAGE PRESSURE RESIDUAL 424084E-04 
STANDARD ERROR OF A SINGLE PRESSURE RESIDUAL 200Z02ZE—03 


# DATA FOR THIS RUN WERE OMITTED FROM THE CALCULATIONS. 


a - : im 
ta - ; 

oe . 
; wi 7 es a. 
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TABLE 9. — PRESSURE RESIDUALS, IN ATM, FOR A QUARTIC FIT IN P 
TO REPRESENT Zy TREATING ALL RUNS SIMULTANEOUSLY 
WEIGHTING FACTOR P SUB R TO THE -3/4 
R=4 


DEV. FROM AVG. 


RUN NO. (P»,OBSe—-P,yCAL. ) PRESSURE RESIDUAL 
| Soe ee ie91L69E-03 2005 15E-03 
HE-0-—2 1.4018E-03 1.5222E-03 
HE-0-3 226712E-04 3-8751E-04 
HE-0-4 403512E-04 5.25551E-04 
HE-O-5 —}..0282E-05 -9.0783E-04 
HE-O-6 —-2.2081E-04 —-1.0042E-04 
HE-6—7 1.2416E-03 1.3620E-03 
HE-0-8 1..795ZE—-03 L23356E-—93 
HE-o-9 —~240332E-04 “Ree JOO Oo 
HE-0O-10 9-1706E-04 1.0374E-03 
HE~O-—L1 leT fi SE~O4 3.1754E-04 
jp, re ee LeiZise-05 1.2422E-03 
RHE -O-¢ 35 facts tE-05 1.4135E-03 
HE-0-14 72.2289E-04 864328E-04 
ipa ames ee ~1.2864E-03 =—1 «1 660E-05 
HE-0-16 -6.6f414E-04 =~525375E-04 
He-O~ty =7 «lL 51l6E-04 —-5.9477E-04 
HE=-0-18 —6 -0621E-04 -4.8582E-04 
HE-O—19 -2.3467E-03 2 et LODE=O3 
HE-O-70 -2e3174E-03 “2 ITOE-O03 
ME=-0-21 —1.2094E-03 -] .0890E-03 
RE-O-22 we he ad ae 1 =1s8T68E-03 


AVERAGE PRESSURE RESIDUAL -1.2039E-04 
STANDARD ERROR OF THE AVERAGE PRESSURE RESIDUAL 208136E-04 
STANDARD ERROR OF A SINGLE PRESSURE RESIDUAL Lew 893-05 


% DATA FOR THIS RUN WERE OMITTED FROM THE CALCULATIONS. 


AY »2JAUQT29R 3RU22 999 - .@ S4BAT 
TAIAT << THB23994R OT . 
au2 4 SOTDAR OMITHOT Sh 


(»jAJ, 9,280 4) 
tO-Jea fe. | 
FOR JB (G8.I 
‘0-4 ‘On-re 
Vg, lees 

45 8S0.!- 


Ties Fans, . 


d4T [ELS 
Pos. f= 


FOE IT f \AtiIdi24i9 3aUe2 Jaq Ae 
fF JAUZEAR® ZJOAR3IVA SHT AO AGAR 
[caA 3ARUZCIAI JjOn1l2 A AD RONAZ CRACRAT 
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TABLE 9. — PRESSURE RESIDUALS, IN ATM, FOR A QUARTIC FIT IN P 
TG REPRESENT Z>y TREATING ALL RUNS SIMULTANEOUSLY, 
WEIGHTING FACTOR P SUB R TO THE -3/4 


R=5 


DEV. FROM AVG. 


RUN NO.» (P,OBSe-P,CAL.) PRESSURE RESIDUAL 
HE=-0-2 9-8296E-04 1.5274E-03 
He-o=2 1.3768E-03 1292P2E-03 
HE-O-3 8.701 7E-04 1.4146E-03 
HE-0-4 9.8514E-04 2.2 296E—-05 
apes) dash -3.8596E-04 1.5850E-04 
HE-O-6 327614E-05 5 28208E-04 
HE~®—7 6 64391E-04 lowlLB83E-03 
HE-O-8 1.3104E-03 1.8548E-03 
as | mg -3.4404E-03 —-2.8960E-03 
HE-0-10 ~1.0876E-05 -5 .4320E-04 
HE 3.01 Z209E-04 1.-f(237E-04 
HE-D-12 —FeZO019E-O4 -1.7571E-04 
+HE-—)-13 8..9225E-04 1.4367E-03 
HE-0-14 205482E-04 7 o9929E-04 
HE-9-i5 -1.4083E-03 —8 «63 90E—-04 
HE~“O—16 -9.0554E-04 -3.6107TE—-04 
HE=)-iv “Be ZOZE<OS =D DT QEAO4 
HE=-0-18 O29 1 3E-O4 -3 ,.8466E-04 
RE=“019 -201431E-03 “4 .DIS6E-03 
Bi) -20 =O ei tT EOS =. b732E-05 
RE=0=21 Haier oes —7o.f845E-04 
HEB 2 —-1.8422E-03 “1 VOTTeE-03 


AVERAGE PRESSURE RESIDUAL 3~-5.4447E-04 
STANDARD ERROR OF THE AVERAGE PRESSURE RESIDUAL 208437E-04 
STANDARD ERROR OF A SINGLE PRESSURE RESIDUAL 1.3031E-03 


# DATA FOR THIS RUN WERE OMITTED FROM THE CALCULATIONS. 
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TABLE 9. — PRESSURE RESIDUALS, IN ATM, FOR A QUARTIC FIT IN P 
TO REPRESENT Z> TREATING ALL RUNS SIMULTANEOUSLY, 
WEIGHTING FACTOR P SUB R TO THE =3/4 


R=6 


DEV. FROM AVG. 


—_—_ 


RUN NO. (P,OBS.—-P,CAL.) PRESSURE RESIDUAL 
1S aaa 5 Lec lse=-63 Pes3Z2GE=03 
Re=0=2 1.2157E-03 E.«Z2899E-03 
HES KeZ599E-65 1.3341E-03 
HE-0O-4 1.4043E-03 1.4785E-03 
HE-@-5 320260E-04 36 1679E=-04 
HE-0-6 6.9078E-04 1 6498E-04 
HE-o=7 220855E-04 208274E-04 
HE—-O-6 Lwi®S7E-OS Pel fove-es 
hE-6+9 -224364E-03 OC SEOs 
HE=-G=10 ~2 2 3606E-04 -~1.6187E-04 
Pre-G—iil 42e6320E-04 503139E-04 
Hib—i—i 2 —2«4841E-04 -1.7421E-04 
e-O-T3 1e5738E-04 8e3157E-04 
HE-0-14 2c01l477E-04 228896E-04 
ee i -503544E-04 —~4.6124E-04 
HE-G—-16 -7T62157E-04 —604737E-04 
Hie-6= 17 -8.0793E-04 =f. gis f FE~Oy 
HE-0-—18 —7 e0474E-04 -~623055E-04 
HEB tsS -1.4243E-03 <heasoOLE=-63 
HE-O= 20 =) s80T2ZE-03 -9.27T02E-04 
HE-O~2 1 -5.1699E-04 -5 .0280E-04 
HE-O- 22 -9.8847E-04 ~9.1428E-04 
AVERAGE PRESSURE RESIDUAL =-7.4194E=-05 
STANDARD ERROR OF THE AVERAGE PRESSURE RESIDUAL 202335E-04 
STANDARD ERROR OF A SINGLE PRESSURE RESIDUAL 1.0235E-03 


# DATA FOR THIS RUN WERE OMITTED FROM THE CALCULATIONS. 


UG A 904 ~MTA WE e@sauaTeaa apuezaae = 22 agen 
1? PIA i140 BAUAITASRI «\ THAZIAVAN oT ; 
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TABLE 9. — PRESSURE RESIDUALS, IN ATM, FOR A QUARTIC FIT IN P 
TO REPRESENT Z>y TREATING ALL RUNS SIMULTANEOUSLY, 
WEIGHTING FACTOR P SUB R TO THE -3/4 


Ae7 
DEV. FROM AVGe 
RUN NO.« (P,0BS.-P,CAL.) PRESSURE RESIDUAL | 
HE-0—1 1.6471E-03 1.2247E-03 
HE-0-2 1.7391E-03 1.3166E-03 
HE-0-3 2.4100E-03 1.9876E=03 
HE-0=4 9.3360E-04 5.1117E-04 | 
HE~0-5 1.7530E-03 1.3305E-03 | 
HE-O0—6 1.7213E-03 1.2988E=03 | 
HE-0-7 6 888 1E-04 2 .663BE-04 ; 
HE~0-8 1.2853E-03 8 »62B7E-04 
HE-0-9 =7,1975E=05 -4.,9440E-04 
HE-0-10 -2.9266E-04 =7.1509E-04 t 
HE-O=-11 2.0025E-04 =2,2217E-04 
HE-0=-12 307806E=-04 =4.4364E-05 f 
#HE-0-13 2.1525E-03 1.7301E=-03 | 
HE-O-14 8.6457E-04 4.4214E=04 
HE-0=15 3.8617E-04 -3.,6258E=-05 
HE-O-16 -4,2284E-04 =8.4528E-04 
HE-O-17 =7 6 5538E=-04 -1.,1778E=-03 
HE-0~18 =2,0436E-03 =2 .4660E-03 
HE-0-19 =2 62245E-04 =6.448B8E-04 
HE-0=20 =-5.7982E=-04 -1.0022E=-03 
HE~O—21 -1.0729E-04 —5 .2972E-04 
HE -0=22 =6 64035E=-04 =1.0627E-03 
AVERAGE PRESSURE RESIDUAL 4,2243E-04 


STANDARD ERROR OF THE AVERAGE PRESSURE RESIDUAL Zo 


STANDARD ERROR GF A SINGLE PRESSURE RESIDUAL 1.0894 


$ DATA FOR THIS RUN WERE OMITTED FROM 


THE CALC! 
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TABLE 10. - Mean pressure residual as_a function 
of expansion number and pressure 


Expansion 
No. 


+0 .935 
PLL 
AAL 
281 
284 
223 
Deke 


I+ t+ 1+ I+ 1+ 


1/ 


=" Mean pressure residual/standard error. 


Average Mean pressure Standard 
pressure residual x 10° error x 103 
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E 


PR yeiots-2-09) "10 -” , atm; "and 


N = 1.994199 + 0.000030. 


Observed experimental pressures, least-squares calculated pres- 
sures that satisfy the Burnett relationship, differences between ex- 
perimental and calculated pressures, relative differences between 
experimental and calculated pressures, the sum of the weighted squares 
of the pressure residuals, and the parameters N, B, C, D, and E are 
contained in table 11. Calculated standard deviations of the param-~ 
eters (designated SN, SB, SC, SD, and SE), the calculated variances 
of the parameters (designated S2N, S2B, S2C, S2D, and S2E), and the 
calculated covariances of the parameters (designated S2BC, S2BD, S2BE, 
S2BN, S2CD, S2CE, S2CN, S2DE, S2DN, and S2EN), are also recorded in 
table ll. 

Variances and covariances of P, BP, CP, DP, and EP for integral 
pressure values are listed in table 12. These values are needed to 
calculate the standard deviation of a calculated compressibility 
factor. 

Table 13 gives calculated values of the compressibility factor 
at integral values of the pressure, along with the calculated standard 
deviation SZ. The results of table 13 show that the calculated stan- 
dard error in the compressibility factor at one atm is less than one 
part in 10°, while at the highest experimental pressure the calculated 
standard error is less than one part in 10°. 

We believe the values given in table 13 are the best values ob- 


tainable from our data. 
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IN ATM, 


VARIANCES AND 
TREATING ALL RUNS SIMULTANEOUSLY, 


REL. 


0.00000E-99 
P.93261E-05 
3269959E=-05 
4.297200E-05 
6.08724E-05 
6-.Z21639E-05 
1.60818E-04 
4e20725E-04 


0.00000E-99 
8 -81056E-06 
2-16994E-05 
4643862E-05 
4.23337E-05 
8.+36339E-05 
1.47908E-04 
4.22738E-04 


0.00000E-99 
2e399T8E-65 
4.2.35926E-06 
3-94176E-O07 
7696980E-06 
5 022294E-05 
1.51459E-04 
5-18790E-04 


0.00000E-99 
fed LOOSE-US 
—ferhhc hoe“ OD 
= Sie 9 IOE—OG 
1229186E-05 
5 -88435E-05 
1.68002E-04 
2e23210E-04 


0.00000E-99 
eet EIT ZE-OS 
—L.6385T21E-O5 
~Lels tose -O5 
—~se0 194 1E-O5 
ego 2 OSE —O5 

3265161E-05 


TABLE 11. - EXPERIMENTAL AND CALCULATED PRESSURES, 
PARAMETERS, STANDARD ERRORS, 
COVARIANCES, 

WEIGHTING FACTOR P SUB R TO THE -3/4 
O DEGREE C HELIUM ISOTHERM 
R Py OBS. Py CAL. P,OBS.-P,CAL. 
RUN NO. HE-O-1 
0 6.4728895E+02 6.4728895E+02 0.00000E-99 
1 2.8125289E+02 2.8124746E+02 5 «43552E-03 
2 1.359684 70E+02 ._1.3167982E+02 4.87179E-03 
S 6-.385353I0E101 6.3850335E+01 3-17479E-03 
4 3.21491659E+01 3.1489742E+01 1.91697E-03 
S i-S661273E+01- 1.-5660290E+01 9-82963E-04 
6 7.48217009E-00 7.8204430E-00 be Z578FE—OS 
7 3.9151145E-00 3.9134673E-00 1.64718E-03 
RUN NO. HE-O-2 
0 6.9039330E+02 .6.9039330E+02 0-00000E-99 
1 2e9760733Et02 2.9760470E+02 2262208E-03 
2 1.23883415E+02 1.3883030E+02 3-84563E-03 
3 6.27203083E+01 6.2£7200100E+01 2e98289E-03 
4 3.23114887E+01 3.3113485E+01 1.40187E-03 
5 1.26462213E+01 1.6460837E+01 1.37680E-03 
6 8.2197040E-00 8.2184882E-00 LeZiS76E-03 
T &elI3S961L1E-00 4.1122220E-00 Ve fS912E-05 
RUN NO. HE-0O-3 ; 
O 7.0Y28236Et02 7.0128236E+02 0-00000E-99 
1 3.0170799E+02 3.0170075E+02 heros Ze Gs 
2 1.24061376E+02 1.4061315E+02 6212973E-04 
S G.8GS3559E+0L 6.8033533E701 20681 72E-05 
— Set SPPSZCOETOL 323517053 E401 2e67126E-04 
5 1.6660572E+01....1.6659701E+01 8.701 72E-04 
6 8.318601L1E-00 8.3173412E—-00 banc 5993E-05 
T 4.21639855E-00 4.1615754E-00 2e41007E-03 
RUN NO. HE-O-4 
0 7.0574520E+02 7.0574520E+02 0.00000E-99 
% Se0S37U3S7TEtOZ2 3.033 7T53STEtO2 —-4.00436E-03 
2 1664133952E+102 1.4134116Et02 -1.64255E-03 
S GeGo (S652E70] -G.8373656E101 -3.99241LE-06 
4 3.23682135E+01 3.23681700E+01 4.235127E-04 
5 1.6741807E+01 1.6740822E+01 9285147E-04 
6 &3590666E-00 .8..3576622E—-00 1.40434E-03 
T 4.1826390E-00 4.1817054E-00 9.33609E-04 
RUN NO. HE-O-5 
Oo G.SrbT9osE+02 6.9/8 7955E+02 0-00000E-99 
1 34.0041472E+02 3.0042229E+02 -7.57023E-03 
2 1.2-4005441E+02 1.4005701E+02 -2.60111E-03 
> Gelfi2ioset+ol Gaff TS630Et01 -1.46191E-03 
S Se5S70G1L91ETOL SeS391Z22ZGEt0O1 -1.02822E-03 
S Pe6S97TSISETOL Le6597699Et01L -3.85967E-04 
6 842868246E-00 8.2865220E-00 3202603E-04 
7 4421479419E-00 4.1461888E-00 1.75300E-03 


4.22620E-04 


ERROR IN P 
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TABLE lL. — EXPERIMENTAL AND CALCULATED PRESSURES, IN ATM; 
PARAMETERS, STANDARD ERRORS; VARIANCES AND 
COVARIANCES, TREATING ALL RUNS SIMULTANEOUSLY; 
WEIGHTING FACTOR P SUB R TO THE -3/4 

O DEGREE C HELIUM ISOTHERM 
R P, OBS. Ps CAL P,OBS.-P,CAL. REL- ERROR IN P 
RUN NO. HE-O-6 
0 7.0298868E+02 7.0298868E+02 0 -00000E-99 0-00000E-99 
1 3.20233845E+02 3.0234131E+02 -2.86047E-03 -9.46117E-06 
2. 1.264088835E+02 1.4089168E+02 -3.33150E-03 ~2236463E-05 
3 648163698E+01 6.8163674E+01 20h Tt ZE—O5 3-47066E-07 
4 3.3579834E+01 3.3580055E+01 ~2.20818E-04 — 6.0 57S 92E-O6 
5 1.6690781E+01 1.6690743E+01 3076142E-05 2 a2 5399E=06 
6 8.3334615E-00 8.3327707E-00 6.90785E-04 8.28929E-05 
7 4.1709998E-00 4.1692785E-00 1. 72130E-—03 4.12683E-04 
RUN NO. HE-O-7 
0 7.0208680E+02 7.0208680E+02 0.00000E-99 0.00000E-99 
1 3.0200812E+02. 3.0200278E+02 5033 FO08E-03 1. 7ETL9E-65 
2 .164074887E+02 144074449E+02 4.3843 7E-03 302¥503E-05 
3 6.68097837E+01 6.2.8094903Et01 2293412E-03 4.30868E-05 
4 3423548005E+01 3.3546763E+01 1.24161E-03 © 3270100E-05 
5 1.6674984E+01 1.6674340E+01 6643915E-04 3086156E-05 
6 8.3248260E-00 8.3246174E-00 2e08550E-04 2.50516E-05 
7 4.1658969E-00 4.1652080E-00 6288819E-04 1.65347E-04 
RUN NO. HE-0O-8 
© {.O061130E+02 7.006T130E+02 0.00000E-99 0.Q00000E-99 
1 3.0144846E+02 3.0144874E+02 ~2.81659E-04 ~9.34355E—-07 
2 1.464050824E+02 1.44050355E+02 4.69212E-03 3an39S9E—-05 
3 6.7985497E4+01 6.7982316E+01 3218078E-03 4.67862E-05 
4 3.3494053E4+01 3.3492257E+01 1679522E-03 5.235984E-05 
5 1.6648795E+01 1.6647484E+01 1.31040E-03 7e87086E-05 
6 8.3723 743E-00 8.3112685E-00 i. 1057 9E-OS 1233029E-04 
7 4.1598289E-00 4.1585436E-00 ln Z2O550E-05 3.08979E-04 
RUN NO. HE-O0O-9 
0 6.8457479E+02 668457479Et+02 0 e00000E-99 0.00000E-99 
1 2.-9540775E402 2.9541010E+02 -2.35287E-03 -7296483E-06 
2 1e63787O047Et02 . 123787382E+02 —-3235022E-03 —2,..4Z99TFE-05 
3 6.6752988E+01 6.6752677E+01 3210920E-04 426577T7E-06 
4 362896560E4+01 3.42896763Et01 -—2.203321E-04 -6.18062E-06 
5 1.6350585E+01 1.6354026E+01 -3.44049E-03 —-2.10420E-04 
6 8.1629533E-00 8.1653898E-00 -—2.43649E-03 —2.0IB482E—-04 
7 4.0856389E-00 4.0857109E-00 ~-7-19757E-05 —1..76167E—-65 
RUN NO. HE-O-10 
O T.O0727559E102 7.0727559E+02 0.00000E-99 0.00000E-99 
1 320394992E+02 3.0394908E+02 8.32446E-04 22¢13876E-06 
2 124159267E+02 1.4159045E+02 Za22l50E—-03 1..56893E-05 
3 6.28490921E+01 6.8490097E+01 8.23304E-04 1.20206E-05 
4 3.3738977E+01 3.3738060E+01 9.17063E-04 Jot LOELE-OS5 
5 1.6767501E+01 1.6768589E+01 -1.08767E-03 -6.48679E-05 
6 8.3712273E-00 8.3714634E-00 -2.36069E-04 -2.82001E-05 
7 4.1883028E-00 4.1885954E-00 -2.92661E-04 ~6.98 7T59E-65 
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TABLE 11. ~- EXPERIMENTAL AND CALCULATED PRESSURES, 
STANDARD ERRORS, 
TREATING ALL RUNS SIMULTANEOUSLY, 
WEIGHTING FACTOR P SUB R TO THE -3/4 


D 


NOWULWMHF OO SJIDNFWNeH O JTIOWNPWN HO NOupwWyrY Oo 


xX YOUR WNH O 


PARAMETERS, 
COVARIANCES, 


P, OBS. 


120340360E+02 
3202498 76E+02 
1.40961 79E+02 
6.8196393E+01 
3 63595563E+01 
1.6697915E+01 
8 .3369834E-00 
421713507E-00 


6 69905324E+02 
3 200866 78E+02 
1.4025209E+02 
6.78653 73E+01 
323435772E+01 
1.6618379E+01 
8.2969111E-00 
421518777E-00 


6 o9853895E+02 
3 2006 756 7E+02 
1.4016953E+02 
6./826865E+01 
323416347E+01 
1.6609979E+01 
8.2927140E-00 
4215003 76E-00 


4 20725396E+02 
320393974E+02 
1.4158487E+02 
6 68487098E+01 
3eS3 T35978BE+O01 
1.6766788E+01 
8.3707330E-00 
4-1888843E-00 


7205 75241E+02 
32033 7578E+02 
1.4134198E+02 
6 .8374213E+01 
3-3681291E+01 
1.6740047E+01 
8.3570058E-00 
4.1813150E-00 


O DEGREE C HELIUM 


Py CAL. 

RUN NO. 
120340360E+02 
320249702E+02 
1.4095937E+02 
6 28195302E+01 
303595366E+01 
1.669828 7E+01 
823365202E-00 
421711504E-00 

RUN NO.» 
6 -9905324E+02 
300086341E+02 
1.4024894E+02 
6./863330E+01 
3¢3434651E+01 
1-6619100E+01 
8.2971595E-00 
401514997E-00 

*RUN NO. 

RUN NO. 
6 o9853895E+02 
30.0067014E+02 
1.4016485E+02 
6¢7824032E+01 
303415624E+01 
1-6609724E+01 
822924993E-00 
42¢1491730E-00 

RUN NO. 
100725396E+02 
320394098E+02 
1.4158693E+02 
6 o8488452E+01 
3e3737264E+01 
1,.6768197E+01 
8.3712684E-00 
4.1884981E-00 

RUN NO. 
120575241E+02 
320337808E+02 
1-4134233E+02 
628374205E+01 
323681966E+01 
1.6740953E+01 
8.3577273E-00 
4.1817379E-00 


ISOTHERM 


P,OBS.—-P,CAL. 


HE-—-O-—-11 


0.00000E-99 
1. 74220E-03 
2241802E-03 
L2O9TL7E-03 
1.97155E-04 
-3.612095E-04 
4063204E-04 
2200252E-04 


HE-O)-i2 


0.O00000E-99 
BedlZ7T9E-03 
3.045862E-03 
2204313E-03 
1.12184E-03 
-7.20190E-04 
-2248410E-04 
3218067E-04 


HE-6-17s 
HE-O-14 


0,.00000E-99 
5253418E-03 
4.67568E-03 
2e83258E-03 
1o22890E-04 
2254821E-04 
201477T4E-04 
8.645 78E-04 


He eS 


0.00000E-99 
=LineoeBOE- O32 
=2a06E25E-03 
~Jo35368E-03 
-~1.28645E-03 
-1.40837E-03 
—-5 235440E-04 
36 86172E~04 


HE-0-26 


0-00000E-99 
~2.29483E-03 
~3.56484E-04 

7266934E-06 
~60174146E-04 
-9.05548E-04 
—~7221573E-04 
-4.22849E-04 
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IN ATM, 


VARIANCES AND 


REL. ERROR IN P 


0.00000E-99 
5» 15936E-06 
totnooae-OS 
1.60004E-05 
5 286849E-06 
—Zel2E59E~O5 
5 255602E-05 
4.80065E-05 


0.00000E-99 
1e212102E-05 
2e25210E-05 
3201056E-05 
sO ed Nees 6 
le 33369E-05 
-2299400E-05 
9..10593E-05 


0.00000E-99 
1.84058E-05 
3<33573E-05 
4e17619E-05 
2o6328E—-05 
1.53414E-05 
2«058992ZE-05 
2008330E-04 


0.00000E-99 
-4,05609E-06 
—-1,.45584E-05 
-1.97654E-05 
“3.8452 9E-05 
—8 -L9STIE-OS 
—6.39657E—05 
9.21898E-05 


0.00000E-99 
-7-256431E-06 
—Zeocea3E-06 
1.12167E-07 
-2.00154E-05 
-5.40947E-05 
-8 .63435E-05 
—-1.01128E-04 


DATA FOR THIS RUN WERE OMITTED FROM THE CALCULATIONS. 
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TABLE 11. = EXPERIMENTAL AND CALCULATED PRESSURES, IN ATM, 
PARAMETERS, STANDARD ERRORS, VARIANCES AND 
COVARIANCES, TREATING ALL RUNS SIMULTANEOUSLY, 
WEIGHTING FACTOR P SUB R TO THE ~3/4 


O DEGREE C HELIUM ISOTHERM 


R P, OBS. P, CAL. P,OBS.-P,CAL. REL- ERROR IN P 
RUN NO. HE-O-17 
OF Ge9911 2336702 ~~ 6.9911233E+02 0.00000E-99 0.COO000E-99 
1 3.0088594E+02 3.0088561E+02 3,.30496E-04 1.09841E-06 
2 1.4025861E+02 1.4025860E+02 le2oi PBE-05 8296330E-08 
3 6.7867489E+01 6.7867845E+01 ~—-3.55787E-04 —-5.24238E-06 
4 3.3436121E+01 3.3436836E+01 -7.15165E-04 —2¢L3890E—-05 
5 1.6619056E+01 1.6620177E+01 —-1.12026E-03 -6.14083E-05 
6 8.2968869E-00 8.2976948E-00 -8.07930E-04 ci Pe te EF Po | asd 8 es. 
7 4.21510115E-00 4.1517669E-00 -7.55380E-04 -1.81974E-04 
RUN NO. HE-O-18 
0 T7.0498797TEt02 7.0498797E+02 0.00000E-99 0-O00000E-99 
1 3.0309166E402 3.0309140E+02 2062479E-04 8.66005E-07 
2 1.24121823E+02 1.4121774E+02 4.90793E-04 3047542E-06 
3 6.8316532E+01 6.8316005E+01 5 226922E-04 7of1296E—-06 
4 323653188E+01 3.3653794E+01 -6.06215E-04 -1.80136E-05 
5 1.6726144E+01 1.6727074E+01 —-9.29137E-04 5.20 DOOE—O5 
6 8.3501239E-00 8.3508286E-00 ~-7/.04/49E-04 —-8.43999E-05 
1 4.1762502E-00 4.1782938E-00 -2.04364E-03 —4.89348E-04 
RUN NO. HE=0=-19 
0 6.9586241E+02 6.9586241E+02 0.00000E-99 0.00000E-99 
1 2.9965869E402 2.9966376E+02 -5.07965E-03 -1.69514E-05 
2 163972241E+02 1.3972691E+02 -4.50444E-03 —3022384E-05 
3 6.7615486E+01 6.7619330E+01 -3.84323E-03 =5 «68395€-05 
4 343314160E+01 3.3316507E+01 ~—-2.34672E-03 —-7 004422E-05 
5 1.6558741E+01 1.6560884E+01 ~-2,.14313E-03 -1,29426E-04 
& S.2Z6679T5E-00 8.2682218E-00 -—1.42434E-03 “Let ce FOE =-04 
7 4.1368299E-00 4.1370523E-00 <~-2.22456E-04 ~5.37746E-05 
RUN NO. HE-0-20 
0 6.8620552E+02 6.,8620552E+02 0.Q00000E=-99 0.00000E-99 
1 2.9602082E+02 2.9602559E4+02 -4.77402E-03 —-1,.61273E-05 
2 123813812E+02 1.3814216E+02 ~4,.04212E-03 -20.92614E-05 
3 6.26875363E401 6.6878222E4+01 ~-2.85917E-03 4.2753 7E-05 
Mw SerS552GZE40L 3Se2957579Et01 —2.231742E-035 =-7.03202E-05 
5 1.6381782E+01 1.6384000E+01 =-2.21771E-03 =P. 353 TTE-O4% 
6 8.1792899E-00 8.1802911E-00 ~-1.00122E-03 =-1.22409E-04 
7 4.20925711E-00 4.0931509E-00 ~-5.79828E-04 -1.41678E-04 
RUN NO. HE-0O-21 
O14 7435393E402  7.27435393E+02 0.00000E-99 0.00000E-99 
Pe eectee OIEtOl 3.288261 8E+02 =-1.211120E-03 “Jed l SSLE-OG 
2 1.25234062E402 1.5234309E+02 ~-2.46583E-03 -1.61863E-05 
3 723499612E+01 7,3499488E+01 1623778E-04 1.68407E-06 
4 3.6158408E+01 3.6159617E+01 -1.20940E-03 = 3.03447 3E-05 
SB £27959522E401 1.7960845E+01 ~-1.32292E-03 ~7.e36615E-05 
6 8.9632908E-00 8.9638678E-00 ~-5.76999E~04 -60.43736E-05 
T 464841902E-00 4.4842975E-00 <-1.07298E-04 “Peo geoleE-O5 
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TABLE 11. —- EXPERIMENTAL AND CALCULATED PRESSURES, 
STANDARD ERRORS, 
TREATING ALL RUNS SIMULTANEOUSLY, 
WEIGHTING FACTOR P SUB R TO THE =—3/4 


y,”) 


NOW PWM O 


PARAMETERS, 
COVARIANCES, 


P, OBS. 


528218590E+02 
225608357E+02 
1.20581 77E+02 
5 .8627782E+01 
208951908E+01 
1.4406838E+01 
16196 7975E—-00 
3 26018337E-00 


O DEGREE C HELIUM 


Py, CAL. 
RUN NO. 
528218590E+02 
225609148E+02 
1.2058508E+02 
5-8630102E+01 
208953845E+01 
1.4408680E+01 
701977859E-00 
3.6024741E-00 


ISOTHERM 


P,OBS.-P,CAL. 
HE-0-22 


0.00000E-99 
7#292628E--03 
=3.S50127E-03 
=Z2ed200Z2E—03 
=L.295TZ1E-03 
-1,.84220E-03 
-9.88478E-04 
-6.40355E-04 
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IN ATMg 


VARIANCES AND 


RELe ERROR IN 

0.-00000E-99 
-3.08894E-05 
=Z2543T19E-05 
=oe959F-20E-05 
S6569L15E-05 
-1.27869E-04 
-1.37349E-04 
—-1,.77785E-04 


SUM OF THE WEIGHTED SQUARES OF THE RESIDUALS 3.54143E-05 


PARAMETERS 
N 1.994198697E-00 
B 52352783000E-04 
C -6.909896 798E-08 
D 20447895123E-11 
E —-7.145536278E-15 
Seo 1 .3139E=15 
SeOR=1.89920E-11 
SZCN -6.07935E—-14 
S2EN 4.90323E-20 


VARIANCES. STANDARD 
S2N 9.01723E~+10 SN 
S$2B 4.52016E-13 SB 
S2C 5236547E-18 Se 
S2D 1.37560E-23 SD 
S2E 4.35240E=-30 Se 
COVARTANCES 
S2BD 2.38348E~-18 S2BE 
SelD —G.52185E=21 s2CeE 
SZ2DE -7268235E-27 S2DN 


ERRORS 


3,00287E-05 
6.72321E-07 
2231634E=-09 
Sed OB9TE-12 
2e08624E=-15 


“Lec ILTSE-21 
4.6893 7E-24 
—9 o21004E-17 


rAuCcE 


HAS ¢ OH 


DAD i A 


+ 
¢) ” 
1) 
rf) 
4) 
fy 
4 - 


7% 


WA 


4Hy O01 


‘TAIZI IAD QUA JAIME A ISR S — otk 
ARF GSAQMAT2 ,2R3TaRiaas 

114A QWMITASAT 
H 68U2 9. FOTIRAS OALTMSTS 


MUPFJIIN 2D S3RosG 


,220Me THA VOD 


SO+sore si 8s-e So+qoceesrTsaae 
Lp . Fg FenaAr _ © oO-+ ~ 8 ad s0a8 rf 


e SQvsd9G28208.) §0+3T) (8d 
fOr 48O1ofes.2 foO+3Sery Sai 

[Or Geehecdes.S 16398001 CR 

. fE2dOS5RGe08.)] Loraseseoaeeee 
OO-3OPSTTCL.T O0-  3apets v 
‘fpf asda. Ho- oT feR Re 


ua 


iC} iAAUIGS GATRO? I aT. 
hE ABN 22373MA RAR Oe 
STig.se Ase. Nd ITO SESE = 
= Chen VO RT 4 

APaE Le PS. yO Deey Ae sege, + 
eeTEut Ake TT" HES PROSTAR 
\OSSB@ELA 3h { «389s acetates 


aL 308 feed Be 
({-20S°°R sae 
ol-32eepo ee 
OS ~3ESS0e.e © 


a1 


TABLE 12. — VARIANCES AND COVARIANCES FOR INTEGRAL PRESSURES, 


PRESSURE, ATM 


1.-000E-00 
2 eO000E-00 
5 -O0O0E-00 
1 O00E+01 
2 e500E+01 
5-O00E+01 
7 e500E+01 
1.000E+02 
LeZ50E+02 
1.500E+02 
2 0O00E+02 
2e500E+02 
3 -0O00E+02 
3-500E+02 
4 .000E+02 
4.500E+02 
5 -O00E+02 
6 eO00E+02 
7 e000E+02 
8 -O000E+02 
9-000E+02 
1-000E+03 


PRESSURE, ATM 


1.000E-00 
2 ~000E-00 
5 eOQO0E-00 
1 .OOOE+01 
22500E+01 
5 .OOOE+01 
72500E+01 
1 .000E+02 
16250E+02 
1 2-500E+02 
2 2000E+02 
2 eD00E+02 
3-000E+02 
3 e500E+02 
4.000E+02 
4 -500E+02 
5 2000E+02 
6 -000E+02 
7e000E+02 
8 .000E+02 
9-000E+02 
1-000E+03 


O DEGREE C HELIUM 


S2P 
3«20062E-09 
Te23117E-69 
1.62005E-08 
2e26992E-08 
4052488E-08 
1.30110E-07 
1.89756E-07 
2253640E-07 
3 296468E-07 
6626075E-07 
1.05897E-06 
1.L2i15&E-06 
1.41314E-06 
3e14773E-06 
62f9395E-06 
iw Eb 1 3GER-85 
1-46683E-05 
1.36301E-05 
4219139E-06 
1e97774E-04 
3212412E-03 
20 35944E-02 


S2DP 
1650047E-16 
2.24 2U83E-1G 
2.59306E-16 
1.68106E-16 

—3.33802E-16 
—7.94516E-16 
—5 .88824E-16 
6.72050E-17 
9.04164E-16 
1.68897E-15 
2.45 798E-15 
1.66480E-15 
—5 .96025E-16 
—3.66120E-15 
-—6 .56706E-15 
—8..32776E-15 
—~8.21527E-15 
—2.41547E-15 
1.01116E-15 
-2.92537E-14 
—1.54144E-13 
-4.71702E-13 


S2BP 
3423666E-11 
4.65290E-11 
5298508E-11 
4.256537E-11 
mows tS LE pt 
—Pwb1 9D TE~EO 
-~1.47700E-10 
-4.33067E-11 

1.04154E-10 

Zetec BOE~-2O 
4.26982E-10 

3234745E-10 
<2 29) SSLE-T1 
-5.64493E-10 
-1.11148E-09 
—1l.49619E-69 
—~Po57800E—-09 
~fel6867E-10 

1.68368E-10 
—326825Z2E—-09 
~2.18574E-08 
~7.00412E-08 


SZEP 
429 2Z294E-20 
“Te ll69SE-19 
—1.35206E~-19 
-8.-47941E-20 
1.83120E-19 
4.18826E-19 
2294408E-19 
-6.82188E-20 
~5~eZ21194E-19 
—-9-23B8136E-19 
~tetZ063E—-18 
=f. S0659E-19 
3294462E-19 
Le SOLSE-18 
See roe ~2s 
4.32074E-18 
4.06615E~-18 
6.42822E-19 
—-5.65463E-19 
1.90494E-17 
9251045E-17 
2284455E-16 


ISOTHERM 


SZCP 
-~1.00403E-13 
-1.42698E-13 
—Psf FETSE-13 
=i.Z111L6E-13 

Berni 5SE—1S 
eto Te—iS 
4e22244E-13 
1.93500E-14 
=Se-rV>O0E—-13 
Min CeVELE—-L2 
=e 56242E—12 
‘Det e suo bt 2 
2263409E~-13 
ery Cat US Sere hd 
41313 1E—-12 
a ke A? ald, Seas MF 
5 e52740E-12 
EelooZ5e—i2 
=66c0294E—15 
1.53547E-11 
8 -60910E-1]1 
Peoss7Te—10 


TREATING ALL RUNS SIMULTANEGUSLY, WEIGHTING FACTOR 
Ov SUS ARETE STHE —3/4 
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TABLE 13. - COMPRESSIBILITY FACTORS AND STANDARD ERRORS 
FOR INTEGRAL PRESSURES, TREATING ALL RUNS 
SIMULTANEOUSLY, WEIGHTING FACTOR P SUB R TO 
THE -—3/4 


O DEGREE C HELIUM ISOTHERM 


PRESSURE, ATM Z SZ 
1.000E-00 1.0005352092E=-00 6.95998E-07 
2 .000E-00 1.0010702803E-00 1.37286E-06 
5 .000E-—00 1.0026746670E-00 3.35316E-06 
1.000E+01 1.0053458975E-00 6.53600E-06 
2-500E+01 1.0133391503E-00 1.53999E-05 
5 .000E+01 1.0265941827E-00 2.82401E-05 
7.500E+01 1.0397672917E-00 3689753E-05 
1 .000E+02 1.0528606047E-00 4.78728E-05 
1.250E+02 1.0658761820E-00 5651702E-05 
1.500E+02 1.0788160172E-00 6.10892E-05 
2.000E+02 1.1044761000E-00 6.96171E-05 
2 6500E+02 1.1298554608E-00 7 «49808E-05 
3.000E+02 12154967635 7E-00 7.84889E-05 
3 .500E+02 1.1798250887E-00 8.12136E-05 
4.000E+02 1.2044392122E-00 8.39787E-05 
4..500E+02 1.2288203267E-00 8.73501E-05 
5 .000E+02 1.2529776809E-00 9.16174E-05 
6 0O0O0E+02 1.3006527435E-00 1.02413E-04 
7.000E+02 1.3475169527E-00 1.12068E-04 
8 .000E+02 1 «3936057118E-00 1.12825E-04 

~ 9,000E+02 1 64389372750E-00 1.11374E-04 
1 .000E+03 1.4835127470E-00 1.84123E-04 
Z= 1+BP+ CP + DP + Dr 
B= (5.352760 -00072) -10 ", atm > 
C = -(6.9099 +0.2316)-10 °, atm ? 
D= (2.4479 +0.3709)-10 7+, atm ® 
E = -(7.15 +2.09)-10 +5, atm * 


Zm 
Hl 


1.994199 +0.000030 
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The results recorded in tables 11 through 13 are given in E for- 
mat. The number to the left of E is to be multiplied by an exponent 
of ten given to the right of E. 

Table 14 contains a comparison of the results obtained when the 
data were fitted simultaneously to a fourth degree polynomial in the 
pressure with the results obtained when the data were treated as sep- 


arate runs. A weighting factor of Kins was used in the calculations. 


COMPARISON OF OUR RESULTS WITH THE WORK OF PREVIOUS INVESTIGATORS 

Values of the second virial coefficient of helium at 0° C, obtain- 
ed by various investigators, are listed in table 15 along with the 
value obtained in this research. 

Compressibility factors obtained from our work were compared as 
a function of pressure with results obtained by other workers. Rela- 
tive percentage deviations from our values are listed in table 16. 

Our compressibility factors are higher than those of Holborn and 
Schultze (28) over the pressure range 0 to 150 atmospheres, but lower 
at 200 atm. The differences are plus 0.04 percent at 75 atm and minus 
0.04 percent at 200 atm. 

The results of this work are consistently higher, 0.3 percent at 
200 atm, than those of Michels and Wouters (35); and are higher, 0.18 
percent at 350 atm and 0.06 percent at 1,000 atm, than those of Wiebe, 
Gaddy, and Heins (51). Also, our results are consistently higher than 
the results published by Miller, Brandt, and Stroud (37), with a maxi- 


mum deviation of 0.048 percent at 200 atm. 
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TABLE 14. - ned from 


Comparison of results obtai 


the separate and simultaneous methods 


Z(L atm) Z(700 atm) N 


Data treated as 
separate runs 1.000541274 +0.000001513 | 1.348033 +0.000211 1.994059 +0.000053 


Data treated 


simultaneously 1.000535209 +.000000696 : 1.34/517 +000 E12 1.994199 +.000030 


Bias .000006065 | -000516 -000140 


Data treated as | 
separate runs 5.413582 +0.085Z10 4 -1.0835 +0.0862 299s OCS ~0). 933; 40 31.50 


Data treated 
simultaneously Jeode2to eG ZI ='.69099 2023.16 24479 +.03709 ~.0725637 6209 


Bias .06 104. Behe Ware, 1.054 ~, 862 


75 
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TABLE 15. - Comparison of the second virial coefficient of helium at 
O° C obtained in this work with the value obtained by 
various other investigators 


Investigator | B, cm? mole ? 
Jacquerod and Scheuer (31-32), 1905-1908 13.45 
Onnes (41), 1908 11.49 
Burt (13), 1910 0.00 
Holborn, Otto, and Schultze (26-28) as given 

by Otto (43), 1915 Lies 
Henning and Heuse (23), 1921 T2550 
Baxter and Starkweather (4), 1926 0.00 
Heuse (24), 1926 8.18 
Boks and Onnes (7-8) as given by Nijhoff (39), 

1928 Lied 
Ellestad and Baxter (20), 1928 biz 
Ellestad and Baxter (21), 1929 9.42 
Heuse and Otto (25), 1929 Led. 
Wiebe, Gaddy, and Heins (51), 1931 11.69 
Keesom and Santen (33), 1933, and Santen (44), 

1934 Liv2z/ 
Burnett (12), 1936 Nit day 
Michels and Wouters (35), 1941 P1237 
Oishi (40), 1942 10.23 
Schneider (45), 1949 EieeZ 
Schneider and Duffie (46), 1949 Lie) 
White, Rubin, Camky, and Johnston (50), 1960 12.08 
Miller, Brandt, and Stroud (37), 1961 11.86 
Canfield, Leland, and Kobayashi (14-16), 1962-1965 12209 
Hoover, Canfield, Kobayashi, and Leland (30), 1964 

[Retreatment of Canfield's (14) original data] 11.96 
Blancett (6), 1966 Di 


This work, 1969 11.997 +0.015 
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TABLE 16. - Comparison of values of the compressibility factor of helium at 0° C obtained in this 


study with values obtained by various other investigators (expressed in terms of a 
percentage relative difference) 


Investigator 1 atm 2 atm 5 atm 
Jacquerod and Scheuer (31-32), 1905-1908 -0.0065 — — 
Onnes (41), 1908 0023 0.0046 0.0116 
Burt (13), 1910 0535 = oe: 
Holborn, Otto, and Schultze (26-28) as given by Otto (43), 1915 0011 0021 0051 
Henning and Heuse (23), 1921 ae ere .0197 ~~ —- 
Baxter and Starkweather (4), 1926 20535 — —_ 
Heuse (24), 1926 .0170 ~~ — 
Boks and Onnes (7-8) as given by Nijhoff (39), 1928 .0012 . 0024 0064 
Ellestad and Baxter (20), 1928 0035 — — 
Ellestad and Baxter (21), 1929 POLS — — 
Heuse and Otto (25), 1929 0011 — — 
Wiebe, Gaddy, and Heins (51), 1931 0014 0028 .0069 
Keesom and Santen (33), 1933, and Santen (44), 1934 0032 0065 .0165 
Burnett G12), 1936 -0011 .0022 0053 
Michels and Wouters (35), 1941 0006 0011 0029 
Oishi (40), 1942 .0079 = = 
Schneider (45), 1949 0008 .0016 0038 
Schneider and Duffie (46), 1949 0010 .0020 0048 
White, Rubin, Camky, and Johnston (50), 1960 -.0004 -.0007 -.0019 
Miller, Brandt, and Stroud (37), 1961 .0006 .0012 .0029 
Canfield (14), 1962 as -.0005 -.0010 -.0025 
Blancett (6), 1966 0004 .0008 .0019 
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TABLE 16. - Comparison of values of the compressibility factor of helium at 0° C obtained in this 
study with values obtained by various other investigators (expressed in terms of a 
percentage relative difference)--Continued 


Investigator 


Jacquerod and Scheuer (31-32), 1905-1908 
Onnes (41), 1908 
gant Aso), 1910 


Holborn, Otto, and Schultze (26-28) as given by Otto (43), 1915 


Henning and Heuse (23), 1921 

Baxter and Starkweather (4), 1926 

Heuse (24), 1926 

Boks and Onnes (7-8) as given by Nijhoff (39), 1928 
Ellestad and Baxter (20), 1928 

Ellestad and Baxter (21), 1929 

Heuse and Otto (25), 1929 

Wiebe, Gaddy, and Heins (51), 1931 

Keesom and Santen (33), 1933, and Santen (44), 1934 
Burnett (12), 1936 

Michels and Wouters (35), 1941 

Oishi (40), 1942 

Schneider (45), 1949 

Schneider and Duffie (46), 1949 

White, Rubin, Camky, and Johnston (50), 1960 
Miller, Brandt, and Stroud (37), 1961 

Canfield (14), 1962 

Blancett (6), 1966 


10 atm 


0.0236 


wOUS9 


25) atm 


0.0611 
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50 atm 
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TABLE 16. - Comparison of values of the compressibility factor of helium at 0° C obtained in this 
study with values obtained by various other investigators (expressed in terms of a 


percentage relative difference)--Continued 


Investigator 75 atm 


Jacquerod and Scheuer (31-32), 1905-1908 ae 
Onnes (41), 1908 — 
Burt (13), 1910 <4 
Holborn, Otto, and Schultze (26-28) as given by Otto (43), 1915 0.0401 
Henning and Heuse (23), 1921 o 
Baxter and Starkweather (4), 1926 — 
Heuse (24), 1926 — 
Boks and Onnes (7-8) as given by Nijhoff (39), 1928 1859 
Ellestad and Baxter (20), 1928 — 
Ellestad and Baxter (21), 1929 — 
Heuse and Otto (25), 1929 — 


Wiebe, Gaddy, and Heins (51), 1931 0847 
Keesom and Santen (33), 1933, and Santen (44), 1934 — 

Burnett (12), 1936 0428 
Michels and Wouters (35), 1941 0640 
Oishi (40), 1942 — 

Schneider (45), 1949 0214 
Schneider and Duffie (46) , 1949 “0357 
White, Rubin, Camky, and Johnston (50), 1960 a 

Miller, Brandt, and Stroud (37), 1961 0325 
Canfield (14), 1962 -.0330 
Blancett (6), 1966 .0210 
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125 atm 


8S 


myeacerr (@)* race 032] “O5at 
erarresd (Te)" ‘1a0s = O73 - OFS) ~* 93350 


WIyyer* prsuge’ sxq ‘2cxo , 33 03% iy 
MPT cS* Bapre’® Cem)” sog josue co (20) Taed 7 : 
Gepsesqe® wy prrire (se) * j / 


ppuoTqsr (F2)" joe) j 6373 eh e 
epT*t0)* Tavs. -_ Z = 


Beeeeye oq Aarcbre (32)° Jay! }er0 0802 risy 
peevere (15)" 123¢ __ = CSE OF0¢e 0218 


mM eug geucey (32)* fag>" seq esuces GF)‘ jar var ss a 
SGeqgek* eve peyue 7)’ 593) Bee ‘Ta23 533 

sock (527° Tasa - “ = 

“SOY PeExpsr ¢$r)’ Tase 7 = 

oq Beers: (34)' rose ™ = 


ee HESS 5S) ee Braet pd wriport (35)* fate . ‘re3a = 
wtih i) * gaye “s . a io 
cette Spee al <a (¢}° Tase . th z= 
HB Wag wemee Sh)" fai : i = 
. es bt = “309g se £t¢ (5e- 59) ge Sracu: pa ocr r> - 
1)" Tat - — =~ 
ees SE PENEX (3I~95)" TAd-1908 


= 


TABLE 16. - Comparison of values of the compressibility factor of helium at 0° C obtained in this 
study with values obtained by various other investigators (expressed in terms of a 


percentage relative difference)--Continued 


Investigator 150 atm 


Jacquerod and Scheuer (31-32), 1905-1908 a 
Onnes (41), 1908 — 
Beret )s + 1910 — 
Holborn, Otto, and Schultze (26-28) as given by Otto (43), 1915 0.0107 
Henning and Heuse (23), 1921 = 
Baxter and Starkweather (4), 1926 — 
Heuse (24), 1926 = 
Boks and Onnes (7-8) as given by Nijhoff (39), 1928 — 
Ellestad and Baxter (20), 1928 — 
Ellestad and Baxter (21), 1929 — 
Heuse and Otto (25), 1929 — 
Wiebe, Gaddy, and Heins (51), 1931 NSE 
Keesom and Santen (33), 1933, and Santen (44), 1934 — 
Burnett ga 1936 a 
Michels and Wouters (35), 1941 1436 
Oishi (40), 1942 
Schneider (45), 1949 — 
Schneider and Duffie (46), 1949 — 
White, Rubin, Camky, and Johnston (50), 1960 — 


Miller, Brandt, and Stroud (37), 1961 0459 
Canfield (14), 1962 -.0597: 
Blancett (6), 1966 0299 


200 atm 


250 atm 
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TABLE 16. - Comparison of values of the compressibility factor of helium at 0° C obtained in this 
study with values obtained by various other investigators (expressed in terms of a 
percentage relative difference)--Continued 


Investigator ; 300 atm 350 atm 400 atm 


Jacquerod and Scheuer (31-32), 1905-1908 = a — 
Onnes (41), 1908 oe = ae 
Burt (13), 1910 ae a a 
Holborn, Otto, and Schultze (26-28) as given by Otto (43), 1915 — — — 
Henning and Heuse (23), 1921 oe oe a 
Baxter and Starkweather (4), 1926 so — ‘ee 
Heuse (24), 1926 a me = 
Boks and Onnes (7-8) as given by Nijhoff (39), 1928 — — — 
Ellestad and Baxter (20), 1928 I a ae 
Ellestad and Baxter (21), 1929 = ned aS 
Heuse and Otto (25), 1929 — 28 id 
Wiebe, Gaddy, and Heins (51), 1931 0.1763 0.1766 iy el) 
Keesom and Santen (33), 1933, and Santen (44), 1934 _ — ~~ 
Burnett (12), 1936 

Michels and Wouters (35), 1941 pee bee aaa 
Oishi (40), 1942 = = =~ 
Schneider (45), 1949 ae = = 
~ Schneider and Duffie (46), 1949 — ae pee 
White, Rubin, Camky, and Johnston (50), 1960 —_ — — 


Miller, Brandt, and Stroud (37), 1961 0417 — — 
Canfield (14), 1962 -.0989 -, 1089 thd 1 
Blancett (6 )4 1966 0281 .0250 0218 
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TABLE 16. - Comparison of values of the compressibility factor of helium at 0° C obtained in this 


study with values obtained by various other investigators (expressed in terms of a 


percentage relative difference)--Continued }}”© 


Investigator 450 atm 500 atm 600 atm 
Jacquerod and Scheuer (31-32), 1905-1908 = — oe 
Onnes (41), 1908 =< ae i 
Bore (13), 1910 = as ee 
Holborn, Otto, and Schultze (26-28) as given by Otto (43), 1915 — -— — 
Henning and Heuse (23), 1921 = — eee 
Baxter and Starkweather (4), 1926 = = a 
Heuse (24), 1926 = set tt 
Boks and Onnes (7-8) as given by Nijhoff (39), 1928 _ — — 
Ellestad and Baxter (20), 1928 es = pa. 
Ellestad and Baxter (21), 1929 — = = 
Heuse and Otto (25), 1929 ibe — =e 
Wiebe, Gaddy, and Heins Col? Ly 3l G16 29 QO. 1569 OE be a 


Keesom and. Santen (33),1933, and Santen (44), 1934 — _ _ 
Burnett (12), 1936 ae = oe 
Michels and Wouters (35), 1941 as — — 
Oishi (40), 1942 < as ie 
Schneider (45), 1949 se sak i 


Schneider and Duffie (46), 1949 — _— — 
White, Rubin, Camky, and Johnston (50), 1960 _ — _ 
Miller, Brandt, and Stroud (37), 1961 — — — 
Canfield (14), 1962 -.1245 -.1303 a 
Blancett. (6), 1966 .0188 .0163 .0134 
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TABLE 16. - Comparison of values of the compressibility factor of helium at 0° C obtained in this 


study with values obtained by various other investigators (expressed in terms of a 


percentage relative difference)--Continued 


Investigator 


Jacquerod and Scheuer (31-32), 1905-1908 

Onnes (41), 1908 

Burt (13), 1910 

Holborn, Otto, and Schultze (26-28) as given by Otto (43), 1915 
Henning and Heuse (23), 1921 

Baxter and Starkweather (4), 1926 

Heuse (24), 1926 

Boks and Onnes (7-8) as given by Nijhoff (39), 1928 
Ellestad and Baxter (20), 1928 

Ellestad and Baxter (21), 1929 

Heuse and Otto (25), 1929 

Wiebe, Gaddy, and Heins (51), 1931 

Keesom and Santen (33), 1933, and Santen (44), 1934 
Burnett (12), 1936 

Michels and Wouters (35), 1941 

Oishi (40), 1942 

Schneider (45), 1949 

Schneider and Duffie (46), 1949 

White, Rubin, Camky, and Johnston (50), 1960 
Miller, Brandt, and Stroud. (37), 1961 

Canfield (14), 1962 

Blancett (6), 1966 


700 atm 


tl 


800 atm 


0.0921 


900 atm 


0.47 37 


1000 atm 


0.0578 
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Our data are consistently lower than Canfield's (14), amounting 
to 0.13 percent at 500 atm, but are only slightly higher than the re- 
sults obtained by Blancett (6). The maximum deviation from Blancett's 
data is 0.031 percent at 200 atm, and the deviation decreases to 


0.014 percent at 700 atm. 


DISCUSSION OF RESULTS 

In all of the various data treatments, we assumed that the pres- 
sure distortion coefficients of the pressure vessels were exactly 
known, that the expansion numbers were free of error, that random 
errors occurred only in the observed pressures, and that these errors 
were normally distributed. We also assumed that all data were taken 
etvexactily.O° C. 

Temperature corrections to the observed pressures, based on the 
observed temperatures reported in table 4 and on the PVT properties 
of helium obtained by Wiebe, Gaddy, and Heins (51), did not change 
any of the calculated quantities by more than seven percent of their 
calculated standard errors. No temperature corrections are applied 
to the data in this report. 

Values for the pressure distortion coefficients reported by 
Briggs and Barieau (10) were used for all data treatments. Completely 
neglecting the correction for distortion of the pressure containers 


would introduce an error of about 0.1 percent into the compressibility 
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factor for helium at 700 atm. The distortion coefficients are known 
with a precision of better than 1 percent; therefore, any error in the 
compressibility factor due to error in the distortion coefficients is 
assumed to be negligible. 

The expansion number is an exact cardinal number only if the 
volume Vz is completely evacuated before expansion from V, .- In”prac- 
tice there was always a residual pressure of about 1.3xl10 © atm or 
less after evacuation. The error due to the residual pressure would 
be about 4 parts in 10° at the lowest experimental pressure. The 
sensitivity of the piston gage was about 50 parts in 10® at the low- 
est pressure; therefore, the errors due to a residual pressure are 
negligible at both high and low pressures. 

National Bureau of Standards calibration of our barometer indi- 
cated that barometric pressures should be accurate to better than 
0.0002 atmosphere or better than 0.01 percent of the lowest experimen- 
tal pressure. 

The piston gage area is accurate to better than 0.01 percent ac- 
cording to the manufacturer. Suppose there is a small systematic er- 
ror in the piston gage area which would lead to a constant percentage 
systematic error in the pressures. The evaluation of a compressibility 
factor is essentially the determination of a PV product divided by the 
PV product of the same mass of gas extrapolated to zero pressure; 
therefore, any constant error of the area of the piston gage cancels 


in the calculation of the compressibility factor. It does lead to an 
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error in the pressure with which we identify the compressibility 
factor. 

An error of 0.01 percent in the pressure would be equivalent to 
errors in the compressibility factor that are between 10 and 25 per- 
cent of the standard errors listed in table 13. We conclude that any 
systematic error in the pressure due to an error in the piston gage 
area is small compared to the random errors in the pressure. 

The maximum change of the piston gage temperature from the cali- 
bration temperature was 5° C; therefore, completely neglecting the 
correction for change of piston gage temperature would introduce er- 
rors of less than 0.01 percent into the pressure measurements. The 
magnitude of the probable error in the piston gage area temperature 
coefficient is not known; however, the error introduced into the pres- 
sure measurements due to a systematic error in the temperature coeffi- 
cient should be less than the error introduced by neglecting the 
correction entirely and should be negligibly small. 

National Bureau of Standards calibration of a similar piston 
gage indicated that the piston gage area pressure coefficient could 
be in error by as much as 30 percent. A 30 percent error in this 
coefficient would introduce a systematic error of about 0.01 percent 
into the compressibility factors at the highest pressures. 

We wanted to determine if our treatment of the experimental data 
had the characteristics of the treatment of data which are known to 


be random. To answer this question, a set of data containing known 
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random errors was constructed. We constructed the random error data 
set by accepting the values of the parameters and of the calculated 
pressures contained in table 11 as the true values. We then applied 
random deviations to these true values by the use of random deviates 
(49). Weighting factors of Ree were assumed. The product of the 
weighting factor and the variance was taken as 2.32 - 10 7’, which 
was obtained from table 8. We thus generated a set of observed pres- 
sures that were known to contain random errors of a size consistent 
with our experimental data. 

These generated data were then treated both as separate runs and 
simultaneously. 

In the simultaneous treatment, the calculated value of B differ- 
ed from the assumed true value by 0.92 times the calculated standard 
error of B. The calculated value of C differed from the true value 
by 1.11 times the calculated standard error of C. The corresponding 
PacelOretoreDawase..2o and for BE, 1.33. The calculated value of N 
differed from the true value by 0.52 times the calculated standard 
error of N. The calculated compressibility factors all differed from 
the true values by about 0.90 times the calculated standard errors. 
The product of the weighting factor and the variance was 2.27 ° 10M. 


which is within 3 percent of the true value. 


The treatment of the generated data as separate runs did not yield 


results that were consistent with the input data for the weighting 


factor. 
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Our experimental data treatment has all of the characteristics 
of the simultaneous treatment of data with known random errors, except 
in two respects. First, there is a significant bias between the re- 
sults obtained when the experimental data were treated as separate 
runs and when the data were treated simultaneously. Second, two out 
of seven of the mean pressure residuals differed from zero by more 
than a standard deviation. 

The stated standard deviations are our best estimates of the 
random errors associated with the quantities recorded in this report. 
The deviations are standard deviations calculated by the usual first 
approximation of the law for the propagation of errors (2, 5, 34). 

We suggest that the standard errors of the compressibility fac- 
tors of table 13 should be multiplied by a factor of two to estimate 


the total random and systematic errors. 
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APPENDIX A 
The compressibility apparatus used in this investigation was 
similar to Burnett's original apparatus. For completeness, the com- 
ponent parts of the apparatus are described. 
The high-pressure jacketed containers used in this investigation 


were designed by Miller (36)2/ | The containers were fabricated to 


1/ 


= See list of references at the end of the main report. 


Miller's specifications by Amarillo Helium Plant machine shop per- 
sonnel. All other major equipment items incorporated into the pres- 
ent apparatus were obtained commercially. 

The high-pressure containers were almost identical in size and 
were surrounded by integral oil-filled pressure jackets. The two 
containers had a volume ratio, (V,+Vz2)/V, , of about 1.99 when connect- 
ed with the various valves and fittings which are component parts of 
V, and Vz. Subsequent to Burnett's original publication, most inves-~- 
tigators (including Burnett) have selected a volume ratio of about 
1.5. The pressure jackets of the containers of this investigation 
were used to determine elastic pressure distortion coefficients as 
described by Briggs and Barieau (10), but the jackets were not used 
in the present study. All of the components constituting volumes V, 


and Vz were assembled in a constant temperature bath; therefore, all 
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of the gas sample confined in volume V, or (V, +V2) was at the bath 
temperature. 

The constant temperature bath used in this investigation was a 


2), 


Chandler Engineering Company’ model 22-1 bath assembly. It was 


= a ee ee a a Fe ee ee 


2/ 


Manufacturers are identified to allow the reader to obtain de- 
tailed descriptive literature of commercially available items. 
This identification should not be construed as Bureau of Mines 


endorsement or recommendation of any particular product or 


manufacturer. 


supplied with a 1,000-watt immersion heater, a 500-watt immersion 
heater controlled by a variable transformer, and a 100-watt immersion 
heater. Heat exchange coils of 1/4-inch outside diameter copper tub- 
ing were built into the bath. Various fluids could be circulated 
through the coils. The bath was supplied with an internal circula- 
ting pump for the bath fluid, and a 1/30 horsepower Aminco electric 
motor-driven propeller type stirrer was added to obtain sufficient 
agitation for good temperature control. Refrigeration was supplied 
to the compressibility bath by pumping fluid through the heat ex- 
change coils from an external tank. The external tank contained a 
750-watt immersion heater, a 500-watt immersion heater controlled by 
a variable transformer, a 100-watt immersion heater, and the refrigera- 


tion coils of a Tecumseh model C2513HTK refrigeration unit. The 
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refrigeration unit was controlled by a relay system actuated by a mer- 
cury thermoregulator. Temperature of the fluid in the external tank 
was controlled within about +0.03° C of the set temperature by the 
on-off action of the refrigeration unit. Fluid was pumped through 

the heat exchange coils of the compressibility bath by a Viking model 
F656G positive displacement pump with a built-in relief valve. There 
was a pressure drop of almost 50 psi through the heat exchange coils. 
The pump discharged about 1 1/2 gallons per minute of fluid and was 
driven by and mounted on a 1/3-horsepower explosion proof electrical 
motor. The pump operated at 1750 revolutions per minute. 

Fine temperature control of the compressibility bath was obtained 
by supplying a small amount of extra refrigeration to the bath and 
controlling the temperature with the 100-watt immersion heater. Power 
to the 100-watt immersion heater was supplied through a Hallikainen 
model 1053A Thermotrol temperature controller with a model 1080A 
temperature-sensing element. The temperature controller was actuated 
by a +0.001° C change in the bath temperature. The Thermotrol tempera-~ 
ture controller can be operated in any one of three modes designated 
by the manufacturer as on-off, proportional, or proportional with re- 
set. For on-off operation, the controller is an alternating current 
Wheatstone bridge with the resistance thermometer sensing element the 
variable arm of the bridge. An adjustable arm of the bridge is used 


to select the temperature setting for bridge balance. A change of 
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resistance of the sensing element produces an alternating current 
voltage which is amplified by a resistance-coupled amplifier. This 
voltage is applied as a bias voltage to a thyratron tube which con- 
trols an enclosed mercury relay, thus turning the heater on and off. 

A negative feedback signal is applied during proportional operation 
to produce a heater-on period which is a function of the bridge unbal- 
ance voltage. A positive feedback circuit with an appropriate time 
constant is added to the negative feedback circuit in order to re- 
store the alternating current bridge unbalance to zero during propor- 
tional with reset operation. The temperature controller proportions 
the heat output by time cycle modulations. The controller was oper- 
ated in the proportional-with-reset mode for all of the experimental 
runs of this investigation. The setting of the temperature controller 
was sensitive to changes in the ambient temperature. It was necessary 
to adjust the setting of the controller occasionally in order to main-~ 
tain a particular bath temperature if there were significant changes 
in the ambient temperature. 

The compressibility bath was filled with a 10 percent by volume 
methyl alcohol in water antifreeze mixture for the 0° C compressibil- 
ity runs. A 50 percent by volume ethylene glycol in water antifreeze 
mixture was circulated from the external tank through the compressibil- 


ity bath heat-exchange coils. 
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Compressibility bath temperatures were measured with a Leeds & 
Northrup Company model 8163, serial No. 1586182 platinum resistance 
thermometer. The thermometer was calibrated by the company in February 
1962 by measurements at the ice point, steam point, boiling point of 
sulphur, and boiling point of oxygen. Constants were supplied by the 
company for use in the Callendar and modified Callendar temperature 
interpolation formulas. Temperatures measured with the platinum re- 
Sistance thermometer were in terms of the International Practical 
Temperature Scale (IPTS). The calibrated thermometer should reproduce 
temperatures on the IPTS within +0.01° C, provided the thermometer is 
used with a calibrated resistance bridge and a reliable value is used 
for the resistance of the thermometer at the ice point. 

A value of 25.54626 40.00004 absolute ohms for the ice point re- 
sistance of platinum resistance thermometer serial No. 1586182 was 
derived from measurements made by Briggs (9) at the triple point of 
water. The measurements were made in this laboratory in June 1964 by 
using a Leeds & Northrup Company model 8069 G-2 Mueller bridge, serial 
No. 1603629. All subsequent resistance measurements for thermometer 
No. 1586182 were made with bridge No. 1603629. A constant current of 
2 milliamperes was supplied through the thermometer for all resistance 
measurements. The bridge was calibrated by the manufacturer in August 
1962. Corrections were provided with sufficient precision to determine 


a resistance greater than 1 ohm to about 2 parts in 10° at the time of 
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the calibration. The bridge can be used to measure resistances over 
the range 0 to 111.111 ohms. The smallest decade step is 0.0001 ohn. 
A temperature change of 0.001° C changed the resistance of the plati- 
num resistance thermometer by about 0.0001 ohm at 0° C. The bridge 
was supplied with a mercury-contact commutator for reversing the con- 
nections to the four-lead platinum resistance thermometer for the 
purpose of minimizing the effect of thermometer lead resistance. 

Critical resistance coils of the bridge were mounted by the manu- 
facturer in a thermally insulated chamber, the temperature of which 
was kept constant by means of an electrical heater controlled by a 
thermoregulator. The manufacturer's literature states that the con- 
Stant temperature chamber holds critical measuring resistors constant 
within 40.01° C. The bridge was calibrated with the built in thermon- 
eter reading 34.5° C. All platinum resistance thermometer measurements 
were made with the bridge thermometer indicating a constant tempera- 
ture of 34.5° C. This arrangement eliminated the need to make temp- 
erature corrections to the bridge resistance readings. 

A Keithley Instruments, Inc., model 149 electronic milli-micro- 
voltmeter was used as a bridge null detector. The milli-microvoltmeter 
had a range of measurements from 0.1 microvolt full scale to 100 milli- 
volts. All thermometer resistance measurements were made when the 
bridge thermoregulator was off. The null detector was rather sensi- 
tive to line voltage disturbances; however, the electronic null de- 


tector eliminated many of the problems associated with the use of a 
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moving-coil galvanometer. 

Gas samples in the pressure containers were separated from oil 
of the pressure-measuring system by a Ruska Instrument Corporation 
model 2416, serial No. 9032 differential pressure cell. The differ- 
ential pressure cell used in this investigation was modified by the 
manufacturer so that it could be immersed in the compressibility bath 
fluid. Oil of the pressure-measuring system was separated from gas 
of the sample system by a thin stainless steel diaphragm located with- 
in the differential pressure cell. The differential pressure cell 
was designed to safely withstand an overpressure of 15,000 psi on 
either side of the diaphragm. The pressure differential across the 
diaphragm was kept as small as possible during operation to insure 
maximum diaphragm life. The position of the diaphragm was established 
by the movement of a stem soldered to the diaphragm. As the diaphragm 
moved, the stem moved a core within the coils of a very sensitive 
differential transformer. Output of the differential transformer was 
indicated by an electronic meter. Maximum sensitivity of the differ- 
ential pressure cell meter was about 0.0001 psi. The meter was zeroed 
before each run with atmospheric pressure applied to both sides of the 
diaphragm. 

The zero of the differential pressure cell changed with changing 
operating pressure. Ruska Instrument Corporation provided calibration 
data to correct for a change of the cell zero due to a change of the 


cell operating pressure. Large overpressures were avoided on either 
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side of the diaphragm in order to prevent a change of the cell zero 
during a compressibility run. It was not possible to apply a numer- 
ical correction for zero shift due to Overpressures; therefore, any 
small overpressures that were necessary during a run were applied con- 
Sistently to the oil side of the diaphragm. The upper chamber of the 
diaphragm cell was oil-filled, and the lower chamber was connected to 
the gas sample volumes. Volumes of the upper and lower chambers of 
the cell were about 29.5 cm® and 0.6 cm®, respectively. Gas pressures 
on the sample side of the diaphragm were balanced on the oil side of 
the diaphragm by an oil-lubricated piston gage. 

The piston gage used in this investigation was a Ruska Instrument 
Corporation model 2400 gage, serial No. 9274, with the high-range 
piston-cylinder assembly. It was calibrated in February 1962 by the 
manufacturer by comparison against Ruska Instrument Corporation master 
gage serial No. 7544. The effective area of the Ruska Instrument Corpora- 
tion master gage was determined at the National Bureau of Standards 
and was reported by the manufacturer to be correct to one partitin uo" 
at 25° C. Since the comparison of gage No. 9274 against master gage 
No. 7544 was carried out with a precision of a few parts per million, 
then the calibration accuracy of gage No. 9274 should be of similar 
accuracy as the master gage. 

Ruska Instrument Corporation supplied calibration data for the 


individual piston gage weights, a value for the effective piston gage 
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area at 25° C and atmospheric pressure, a coefficient to correct for 

a change of effective piston gage area as a function of pressure, and 
a coefficient to correct for change of gage temperature from the cali- 
bration temperature. The range of the piston gage used in this inves- 
tigation was from 30 to 12,140 psig. A special weight set was purchased 
with the gage so that the entire pressure range could be covered in 
increments of 0.001 psig. The minimum resolution of the piston gage 
was stated by the manufacturer to be 5 parts per million at full load, 
decreasing to 50 parts per million at tare load. The cylinder of this 
piston gage is of the re-entrant type in which the effective portion 
of the cylinder is subjected to the measured pressure at the external 
cylinder wall; therefore, the gage can be used at high pressures with- 
out an excessive fall rate of the piston and weight load. 

The height of the piston gage was adjusted, by placing accurately 
machined brass blocks under the piston gage legs, so that the refer- 
ence plane of the piston gage was at the level of the diaphragm of 
the differential pressure cell. The zero of a Ruska Instrument Cor- 
poration model 2409 oil manometer was set at the level of the refer- 
ence plane of the piston gage and the diaphragm of the differential 
pressure cell. No oil pressure head corrections were required for the 
apparatus used in this investigation. 

The pressure containers were placed horizontally in the compres- 


sibility bath. The center of the jacketed container portion of V, was 
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about 7.0 cm below the diaphragm of the differential pressure cell, 
and the center of the jacketed container portion of Vg was about 10.2 
cm below the level of the diaphragm. A correction would obviously 
have to be made for a gas head if one were interested in the pressure 
at some level in the pressure containers below the level of the 
diaphragm. The gas head correction is small for helium (less than 1.5 
parts per million at the center of Vg at 800 atmospheres) but would 
be more significant for a gas such as nitrogen. No gas head correc- 
tions were applied to the data of this report. 

Barometric pressures were measured with a Henry J. Green Instru- 
ments, Inc., model 16, serial No. 13346 Fortin type mercurial barometer. 
Barometer serial No. 13346 was calibrated at the National Bureau of 


Standards in May 19642! All barometer readings for the data of this 


eer rene ene re en 
3/ 


National Bureau of Standards Report of Calibration Identification 


No. P7068. On file in the Bureau of Mines Helium Research Center 


library. 


report were made to the nearest 0.05 millimeter of mercury. The 
barometer was supplied with an integral mercury-in-glass thermometer 
graduated in 0.05° C increments. 

Relative humidity in the PVT laboratory was measured with a sen- 


sitive hygrometer manufacturer by Bacharach Industrial Instrument 
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Company. The hygrometer was read to the nearest 1 percent. The 
manufacturer stated an accuracy of +1.5 percent relative humidity for 
the hygrometer. 

The PVT apparatus was evacuated with a Welch Scientific Company 
model 1400B high-vacuum pump. The manufacturer claimed an ultimate 
vacuum capability of 0.1 micron for the vacuum pump. 

Vacuum was measured with a Consolidated Vacuum Corporation model 
GPH-100A vacuum gage with a GPH-001 cold cathode ionization tube. 

The vacuum gage was designed to measure pressures over the range 1xl0~” 
to 25x10 * millimeters of mercury. The actual vacuum attained in the 
compressibility apparatus, as indicated by the vacuum gage, ranged 
between 5 and 10 microns. 

The compressibility containers were filled with gas by use of a 
onivtn model B2C1000 diaphragm-type compressor. The compressor was 
designed for use to 15,000 psi with a compression ratio of 20 to one. 
Less than 5 minutes were required to fill the pressure containers 


to the maximum operating pressure with the Corblin compressor. 
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APPENDIX B 

The following description outlines the procedure used to obtain 
the data. 

Temperature of the compressibility bath was lowered to 0° C by 
pumping refrigerated fluid through the heat exchange coils. The bath 
temperature controller and the refrigerator mercury thermoregulator 
were adjusted until the compressibility bath temperature was maintain- 
ed at the desired level. Temperature of the external tank was con- 
trolled at about -3.5° C in order to maintain the bath at 0° C. 

The gas side of the differential pressure cell was opened to at- 
mospheric pressure. The valve to the oil manometer was opened, and 
the level of the oil in the manometer was adjusted to correspond to 
the level of the diaphragm of the differential pressure cell. The 
readout meter of the differential pressure cell was set at maximum 
sensitivity, and the meter was zeroed. 

The gas system was evacuated, then purged and filled to supply 
cylinder pressure with test gas. Volume V, , which had been previous- 
ly evacuated, was filled with test gas by slowly expanding the gas 
from Vz into V, while simultaneously maintaining a slight overpressure 
on the oil side of the differential pressure cell diaphragm. The 
valve between V, and Vg was closed, then Vz was filled with test gas 
to a high pressure by use of the compressor. Volume V, was filled 
with test gas at high pressure by expansion from Vz while simulta- 


neously maintaining a small overpressure on the oil side of the 
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differential pressure cell diaphragm by use of oil displacement pump 
No. 1. The valve was closed between V, and Vg, and gas in Vz was 
vented to the atmosphere. Volume Vz was then evacuated with the vac- 
uum pump. The rather awkward and inefficient filling procedure de- 
scribed above was necessary because the compressor and compressor 
controls were located in one room and the remainder of the apparatus 
was located in another room. In addition, it was necessary to avoid 
large overpressure on either the oil or gas sides of the differential 
pressure cell diaphragm to prevent a change in the cell zero during 

a run. 

Gas in V, was allowed to reach temperature equilibrium. Gas 
pressure in V, was balanced with oil pressure by adjusting the piston 
gage load until the differential pressure cell meter indicated the 
null position. Temperature equilibrium was assumed to be established 
when the piston gage indicated a stable pressure. The designations of 
weights on the piston gage, piston gage temperature, barometer scale 
reading, barometer temperature, relative humidity, and the resistance 
of the platinum resistance thermometer were recorded. The differen- 
tial pressure cell was closed off from the oil pump and piston gage, 
and the piston gage was lowered to its rest position by slowly de- 
creasing the pressure with the oil displacement pump. The piston gage 
was isolated from the system by closing the valve between the oil pump 


and the piston gage. Oil pressure was raised to equal the gas pressure 
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by use of the oil displacement pump, and the valve to the differen- 
tial pressure cell was opened. Vacuum in Vag was measured with the 
vacuum gage, the valve to the vacuum pump was closed, and the high- 
pressure exhaust valve from Vz was closed. Gas in V, was slowly ex- 
panded into the evacuated Vz. A small overpressure was maintained on 
the oil side of the differential pressure cell diaphragm during the 
expansion. It was important that any over-pressure be maintained 
consistently on the oil side of the diaphragm. If large overpressures 
had been applied alternately to the oil and gas sides of the diaphragm, 
there would have been a change in the zero position of the diaphragm. 
Zero shift of the diaphragm due to large overpressures, alternately on 
the oil and gas sides of the diaphragm, could have introduced errors 
of as much as 0.02 psi in the measured pressures, according to the 
manufacturer's calibration report. 

Gas in V, and Vz was allowed to reach thermal equilibrium, as 
determined by a stable pressure indication. Usually, temperature 
equilibrium was attained in 30 minutes or less. The piston gage was 
isolated from the system, and the differential pressure indicator was 
maintained in the zero position with the oil displacement pump as the 
expansion valve between V, and Vz was slowly closed. Gas in Vo was 
vented to the atmosphere and Vz was evacuated after the expansion 
valve was closed. There was no detectable change of temperature of 


V, or of the compressibility bath temperature during the venting and 
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evacuation of V3. The valve to the piston gage was reopened. Weights 
on the piston gage were adjusted until a null position was indicated 
by the differential pressure cell meter. All of the necessary obser- 
vations were recorded. Then another expansion was made from V, into 
an evacuated Vz. The measuring and expansion sequence was continued 
until the lowest pressure that could be measured with the piston gage 


was reached. 
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